Glueballs at Finite Temperature in SU{3) Yang-Mills Theory 
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Thermal properties of glueballs in SU (3) Yang- Mills theory are investigated in a large temperature 
range from 0.3Tc to l.QTc on anisotropic lattices. The glueball operators are optimized for the 
projection of the ground states by the variational method with a smearing scheme. Their thermal 
correlators are calculated in all 20 symmetry channels. It is found in all channels that the pole 
masses Mq of glueballs remain almost constant when the temperature is approaching the critical 
temperature from below, and start to reduce gradually with the temperature going above Tc- The 
correlators in the 0^"*", and 2^"*" channels are also analyzed based on the Breit-Wigner Ansatz 
by assuming a thermal width F to the pole mass ojq of each thermal glueball ground state. While 
the values of ujo are insensitive to T in the whole temperature range, the thermal widths F exhibit 
distinct behaviors at temperatures below and above Tc. The widths are very small (approximately 
few percent of ujq or even smaller) when T < Tc, but grow abruptly when T > Tc and reach values 
of roughly F ~ a;o/2 at T 1.9Tc. 

PACS numbers: 12.38.Gc, 11. 15. Ha, 14.40.Rt, 25.75.Nq 



I. INTRODUCTION 

The past two or three decades witnessed intensive and 
extensive studies on the t)hase transition of quantum 
chroniodynamics(QCD) [J], which is believed to be the 
fundamental theory of strong interaction. Based on the 
two characteristics of QCD, namely the conjectured color 
confinement at low energies and the asymptotic freedom 
of gluons and quarks at high energies, QCD at finite tem- 
perature is usually described by two extreme pictures. 
One is with the weakly interacting meson gas in the 
low temperature regime and another is with perturba- 
tive quark gluon plasma (QGP) in the high tempera- 
ture regime. The two regimes are bridged by a decon- 
flnement phase transition (or crossover). The study of 
the equation of state shows that the perturbative pic- 
ture of QGP can only be achieved at very high temper- 
atures T > 2Tc. In other words, the dynamical degrees 
of freedom up to the temperature of a few times of Tc 
are not just the quasifree gluons and quarks 0]. Some 
other theoretical studies also support this scenario and 
conjecture that in the intermediate temperature range 
above there may exist different types of excitations 
corresponding to different distance scales [1, |3| render- 
ing the thermal states much more complicated. Apart 
from the quasifree quarks and gluons at the small dis- 
tance scale, the large scale excitations can be effective 
low-energy modes in the mesonic channels as a result 
of the strongly interacting partons Q . The properties of 
the interaction among quarks and gluons at low and high 
temperatures can be studied with thermal correlators. 



There have been many works on the correlators of char- 
monia at finite temperature. Phenomenological studies 
predicted the binding between quarks is reduced to dis- 
solve J/tp at temperatures close to Tc and proposed the 
suppression of charmonia as a signal of QGP [6, 7]. For 
example, potential model studies show that excited states 
like Tp and Xc are dissociated at Tc, while the ground 
state charmonia J/^ and r]c survive up to T = l.lTc [1- 
[isj . However, it is unclear whether the potential model 
works well at finite temperatures p^ . In contrast, many 
recent numerical studies indicate that J/^* and rjc might 
still survive above 1.5Tc p"5l - [l9j . Of course, it is pos- 
sible that the cc states observed in lattice QCD are 
just scattering states. A further lattice study on spa- 
tial boundary-condition dependence of the energy of low- 
lying cc system concludes that they are spatially local- 
ized (quasi)bound states in the temperature region of 
1.11 ~ 2.07Tc [l^]. Obviously, the results of numerical 
lattice QCD studies are coincident to the picture of the 
QCD transition in the intermediate temperature regime. 

Until now most of the lattice studies on hadronic corre- 
lators are in the quenched approximation. Because of the 
lack of dynamical quarks in quenched QCD the binding 
of quark-antiquark systems must be totally attributed 
to the nonperturbative properties of gluons, which are 
the unique dynamical degree of freedom in the theory. 
Since glueballs are the bound states of gluons, a natural 
question is how glueballs respond to the varying tem- 
peratures. At low temperature T ^ 0, the existence 
of quenched glueballs have been verified by extensive 
lattice numerical studies, and their spectrum are also 
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established quite well |21| - |28l |. An investigation of the 
evolution of glueballs versus the increasing temperature 
is important to understand the QCD transition [2^ [s^l 
and the hadronization of quark-gluon plasma (3l| . From 
the point of view of QCD sum rules, glueball masses are 
closely related to the gluon condensate. Lattice stud- 
ies [12] and model calculations indicate that the 
gluon condensate keeps almost constant below Tc and 
reduces gradually with the increasing temperature above 
Tc- Based on this picture, it is expected intuitively that 
glueball masses should show a similar behavior also until 
they melt into gluons [S^l- In fact, there has already been 
a lattice study on the scalar and tensor glueball proper- 
ties at finite temperature [ssj . In contrast to the expecta- 
tion and the finite T behavior of charmonium spectrum, 
it is interestingly observed that the pole-mass reduction 
starts even below Tc (mG(T ~ Tc) ~ 0.8tog(T - 0)). It 
is known that the spatial symmetry group on the lattice 
is the 24-element cubic point group O, whose irreducible 
representations are R = Ai, A2, E, Ti, and T2. Along 
with the parity P and charge conjugate transformation 
C, all the possible quantum numbers that glueballs can 
catch are R^'^ with PC = ++, — h, H — ,, and -I— 1-, which 
add up to 20 symmetry channels. Motivated by the dif- 
ferent temperature behaviors of cc systems with different 
quantum numbers, we would like to investigate the tem- 
perature dependence of glueballs in this paper. 

Our numerical study in this work is carried out on 
anisotropic lattices with much finer lattice in the tem- 
poral direction than in spatial ones. In order to explore 
the temperature evolution of glueball spectrum, the tem- 
perature range studied here extends from 0.3Tc to 1.9Tc, 
which is realized by varying the temporal extension of 
the lattice. Using anisotropic lattices, the lattice param- 
eters are carefully determined so that there are enough 
time slices for a reliable data analysis even at the high- 
est temperature. In the present study, we are only inter- 
ested in the ground state in each symmetry channel R^'~^ . 
For the study optimized glueball operators that couple 
mostly to the ground states are desired. Practically, these 
optimized operators are built up by the combination of 
smearing schemes and the variational method [2l| - |23| . 
In the data processing, the correlators of these optimized 
operators are analyzed through two approaches. First, 
the thermal masses Mq of glueballs are extracted in all 
the channels and all over the temperature range by fit- 
ting the correlators with a single-cosh function form, as is 
done in the standard hadron mass measurements. Thus 
the T-cvolution of the thermal glueball spectrums is ob- 
tained. Secondly, with respect that the finite tempera- 
ture effects may result in mass shifts and thermal widths 
of glueballs, we also analyze the correlators in A'^'^ , A^^ , 
E~^~^, and T2~^ channels with the Breit-Wigner Ansatz 
which assumes these glueball thermal widths, say, change 
Mo into LjQ — iT in the spectral function (see below). As 
a result, the temperature dependence of wq and T can 
shed some light on the scenario of the QCD transition. 

This paper is organized as follows. In Sec. II, a descrip- 



tion of the determination of working parameters, such as 
the critical temperature Tc, temperature range, and lat- 
tice spacing a^, as well as a brief introduction to the vari- 
ational method is given. In Sec. Ill, after a discussion of 
its feasibility, the results of the single-cosh fit to the ther- 
mal correlators are described in details. The procedure of 
the Breit-Wigner fit is also given in this section. Section 
IV gives the conclusion and some further discussions. 



II. NUMERICAL DETAILS 

For heavy particles such as charmonia and glueballs, 
the implementation of anisotropic lattices is found to be 
very efficient in the previous numerical lattice QCD stud- 
ies both at low and finite temperatures. On the other 
hand, the Symanzik improvement and tadpole improve- 
ment schemes of the gauge action are verified to have bet- 
ter continuum extrapolation behaviors for many physical 
quantities. In other words, the finite lattice spacing ar- 
tifacts are substantially reduced by these improvements. 
With these facts, we adopt the following improved gauge 
action which has been extensively used in the study of 
glueballs 
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where /? is related to the bare QCD coupling constant, 
^ = as/at is the aspect ratio for anisotropy (we take ^ = 5 
in this work), Us and ut are the tadpole improvement 
parameters of spatial and temporal gauge links, respec- 
tively, flc = J2c |-ReTr(l - Wc), with Wc denoting 
the path-ordered product of link variables along a closed 
contour C on the lattice, flgp includes the sum over all 
spatial plaquettes on the lattice, Qtp includes the tem- 
poral plaquettes , Vlsr denotes the product of link vari- 
ables about planar 2x1 spatial rectangular loops, and 
ilstr refers to the short temporal rectangles (one tempo- 
ral link, two spatial). Practically, Ut is set to 1, and Us is 
defined by the expectation value of the spatial plaquette. 



A. Determination of critical temperature 

Since the temperature T on the lattice is defined by 
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where Nt is the temporal lattice size, T can be changed 
by varying either Nt or the coupling constant (3 which 
is related directly to the lattice spacing. In order for 
the critical temperature to be determined with enough 
precision, for a given Nt — 24, we first determine the 
critical coupling f3c, because f3 can be changed continu- 
ously. The order parameter is chosen as the susceptibility 
Xp of Polyakov line, which is defined as 
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TABLE I: The simulation parameters for the determination 
of the critical point. The configurations are selected every ten 
sweeps. 

f} Total configurations Thermalization Bin size 



2.80 
2.805 
2.81 
2.815 
2.82 



20000 
30000 
30000 
20000 
8000 



5000 
10000 
10000 
5000 
3000 



1000 
1000 
1000 
1000 
500 




2.S00 2,802 2,804 2,806 2,808 2,810 2,812 2,814 2,816 2,8 



FIG. 1: The xp extrapolation based on the spectral density 
method. The open triangles denote the simulated values of 
Xp, while the filled squares are the extrapolated values. The 
peak position gives the critical /3c ~ 2.808. 



where Q is the Z{3) rotated Polyakov line, 

RePexp[-27ri/3]; arg P e [tt/B, tt) 
e = i ReP; argPe [-7r/3,7r/3) ,(4) 

RePexp[27rV3]; arg P e [-tt, -7r/3) 

and P represents the trace of the spatially averaged 
Polyakov line on each gauge configuration. 

After a /3-scanning on = 24"* anisotropic lattices 
with ^ 5, the critical point is trapped in a very nar- 
row window G [2.800, 2.820]. In order to determine T^, 
more precisely, a more refined study is carried out in the 
/3 window mentioned above with much larger statistics 
through the spectral density method. Practically, the 
spectral density method [3B,|37| is applied to extrapolate 
the simulated xp's at /3 = 2.805, 2.810, and 2.815. In ta- 
bleUare the numbers of heat-bath sweeps for each f3. The 
extrapolation results are illustrated in Fig. [1] where the 
open triangles denote the simulated values of xpi while 
the filled squares are the extrapolated values. Finally, 
the peak position gives the critical coupling constant 
(3c = 2.808, which corresponds to the critical tempera- 
ture Tc « 0.724rn ^ = 296 MeV with the lattice spacing 
ro/a, = 3.476 HI and r^^ = 410(20) MeV. 

With Tc fairly determined, the working coupling con- 
stant /? is set based on two requirements. First, the spa- 
tial volume of the lattice should be large enough in order 
for the glueballs to be free of any sizable finite volume 
effects. Secondly, we require that temporal lattice has a 



TABLE II: Listed are the parameters used to check the critical 
behavior for /3=3.2. The configurations are selected every ten 
sweeps. 

Nt Total configurations Thermalization < P > xp 



60 
48 
44 
40 
36 
32 
28 
24 
20 



2000 
2000 
8000 
8000 
8000 
3000 
2000 
2000 
2000 



500 


-8.73x10- 


-5 


6.65x10 


500 


6.01x10" 


5 


1.81x10 


2000 


2.25x10" 


3 


3.12x10 


2000 


1.72x10" 


2 


9.14x10 


2000 


5.21x10" 


2 


3.10x10 


1000 


8.51x10" 


2 


2.23x10 


500 


0.1253 




2.00x10 


500 


0.1817 




2.09x10 


500 


0.2571 




1.82x10 



0.014 
0.012 
0.01 
0.008 
0.006 
0.004 
0.002 
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FIG. 2; Xp is plotted versus Nt a.t (3 — 3.2. There is a peak 
of Xp near Nt — 40. 



good resolution even at the temperature T ^ 2Tc. Prac- 
tically the working coupling constant is finally set to be 
/3 = 3.2. The lattice spacing at this (3 is set by calcu- 
lating the static potential V{r) on an anisotropic lattice 
24^^ x 128. With the conventional parametrization of V(r), 



V{r) = Vo + ar + 



(5) 



the lattice spacing Og is determined in the units of ro to 
be 



1.65- 



0.1825(7) 



(6) 



where ro is the hadronic scale parameter. If we take 
= 410(20) Afey, we have = 0.0878(4) fm. The 
spatial volume at L = 24 is therefore estimated to be 
(2.1fm)3. On the other hand, using = 296MeV ob- 
tained at (3 = 2.808 as a rough estimate of and ignor- 
ing the systematic error due to finite lattice spacings, Tc 
and 2Tc at (3 = 3.2 are expected to be achieved around 
Nt ^ 40 and Nt ~ 20, respectively. Obviously, the above 
two requirements are all satisfied. 

Based on the discussions above, with a fixed (3 = 3.2, 
the calculations of the thermal correlators of glueballs are 
carried out on a series of lattice 24^ x Nt with Nt = 20, 



4 



r 




FIG. 3: Prototype Wilson loops used in making the smeared 
glueball operators [2^. 



24, 28, 32, 36, 40, 44, 48, 60, 80, and 128, which cover the 
temperature range 0.3Tc < T < 2Tc. As a cross-check, 
Xp at different Nt are calculated first and the results are 
shown in Fig. [2] and Table |lTl It is clear that the ex- 
pectation value of the Polyakov line drops to zero near 
Nt = 40 and the peak position of xP: which gives the 
critical temperature, is trapped between Nt — 36 and 
Nt — AO. In practice, we do not carry out a precise de- 
termination of Tc at /3 = 3.2, but take the temperature 
at Nt = 38, T w (38at)-i = (38a^/0"^ = 296MeV, as 
an approximation of Tc (/3 — 3.2), to scale the temper- 
atures involved in this work. It should be noted that, 
owing to the lattice artifact, the critical temperature Tc 



determined at different lattice spacing (or (3) may differ 
from each other. The closeness of Tc(/3 = 2.808) and 
Tc{P = 3.2) may signal that the lattice spacing depen- 
dence of Tc is mild in this work due to the application of 
the improved gauge action. 



B. Variational method 

It is known that many states contribute to a hadronic 
two-point function. Ideally one can extract the informa- 
tion of the lowest-lying states from the two-point function 
in the large time region if it lasts long enough in the time 
direction. This is the case for some light hadron states, 
such as TT meson, K meson, etc. However, for heavy par- 
ticles, especially for glueballs whose correlation function 
are much more noisy than that of conventional hadrons 
made up of quarks, their two-point functions damp so 
fast with time that they are always undermined by noise 
rapidly before the ground states dominate. Practically, 
in the study of the glueball sector, in order to enhance the 
overlap of the glueball operators to the ground state, the 
commonly used techniques are the smearing schemes and 
the variational techniques. In this work, we adopt the 
sophisticated strategy implemented by the studies of the 
zero-temperature glueball spectrum [2l| - |23j , which is out- 
lined below. First, for each gauge configuration, we per- 
form six smearing/fuzzing schemes to the spatial links, 
which are various combinations of the single-link proce- 
dure (smearing) and the double-link procedure (fuzzing) 



U-^ix) = Psu(3){Ujix) + Xs Uk{x)U,ix + k)ul{x + j)}, 

±(Mj) 

Ujix) = Psu{3){UAx)U,{x+j) + Xf Uk{x)U, {x + k)U,{x+ J + k)Uk{x + 2j)}, (7) 
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where Psu(3) denotes the projection into SU{3) and is 
realized by the Jacobi method [s^. The six schemes are 
given explicitly as s\°, s\^^, sf^, hf<S)si°^ hf<S)s\l, 
f\f ^ where s/ f denotes the smearing/fuzzing pro- 
cedure defined in Eq. (O, and Xs/Xf the tunable pa- 
rameter which we take As = 0.1 and A/ — 0.5 in this 
work. Secondly, we choose the same prototype Wilson 
loops as that in Ref. [1^ (as shown in Fig. such 
that for each smearing/fuzzing scheme, all the different 
spatially oriented copies of these prototypes are calcu- 
lated from the smeared gauge configurations. Thus for 
a given irreducible representation R of the spatial sym- 
metry group O, say, R = Ai, A2, E,Ti, or T2, a real- 
ization of R can be a specific combination of differently 
oriented Wilson loops generated from the same proto- 
type loop (one can refer to Ref. [23l| for the concrete 



combinational coefficients). The glueball operators (j) 
with the quantum number R^'~^ are thereby constructed 
along with the spatial refiection and the time inversion 
operations. In practice, we establish four realizations 
of each R^'~^ which are based on four different proto- 
types, respectively. Therefore, along with the six smear- 
ing/fuzzing schemes, an operator set of the same specific 
quantum number i?^*^ is composed of 24 different op- 
erators, {0Q, a — 1,2,..., 24}. The last step is the im- 
plementation of the variational method (VM). The main 
goal of VM is to find an optimal combination of the set of 
operators, •l' = ^a'/'ai which overlaps most to a specific 
state (in this work, we only focus on the ground states). 
The combinational coefficients v = {va, a = 1, 2, . . . , n} 
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can be obtained by minimizing the effective mass, 

m{tD) = -—\n^ . , (8) 

at ti) — I, where C'afj (t) is the correlation matrix of the 
operator set, 



^a^W = I](O|0a(i + r)<^^(r)|O). 



(9) 



This is equivalent to solving the generalized eigenvalue 
equation 



(10) 



and the eigenvector v gives the desired combinational co- 
efficients. Thus, the optimal operator that couples most 
to a specific states (the ground state in this work) can be 
built up as 



$ = ^ Va(j)a, 



(11) 



whose correlator C{t) is expected to be dominated by the 
contribution of this state. 



III. DATA ANALYSIS OF THE THERMAL 
CORRELATORS OF GLUEBALLS 

AU 20 R^^' channels, with R = Ai, A2, E,Ti,T2 

and PC = ++,H — +, , are considered in the 

calculation of the thermal correlators of glueballs on 
anisotropic lattices mentioned in Sec. II. At each tem- 
perature, after 10000 pseudo-heat-bath sweeps of ther- 
malization, the measurements are carried out every three 
compound sweeps, with each compound sweep composed 
of one pseudo-heat-bath and five micro-canonical over- 
relaxation(OR) sweeps. In order to reduce the possible 
autocorrelations, the measured data are divided into bins 
of the size rimb = 400, and each bin is regarded as an in- 
dependent measurement in the data analysis procedure. 
The numbers of bins A^bin and n,„b at various tempera- 
ture are listed in Table Hill 

Theoretically, under the periodic boundary condi- 
tion in the temporal direction, the temporal correlators 
C{t, T) at the temperature T can be written in the spec- 
tral representation as 



^ |(n|$|m)|2 f En 
= ^TTTTTTT^f^ exp 



2Z{T) 



En 



2T 



X cosh 



dujp{uj)K(uj, T), 



(12) 
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400 
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400 
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0.86 


400 
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36 
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400 


40 


32 


1.19 


400 


56 


28 


1.36 


400 


40 


24 


1.58 


400 


40 


20 


1.90 


400 


40 



with a T-dcpendent kernel 

cosh(w/(2T) - wt) 



sinh(w/(2r)) 



and the spectral function, 
^(-) = L 2Z{T) ' 

X (,5(W - {En - Era) - Siuj - {En, - S„)),(14) 

where Z{T) is the partition function at T, and £"„ the 
energy of the thermal state \n) (|0) represents the vac- 
uum state). In the zero-temperature limit (T — 0), due 
to the factor exp(— £^„i/r), the spectral function p{lli) 
degenerates to 



p(^) = E' 



2Z{0) 



{S{uj - En) - S{u; + En)) , (15) 



thus we have the function form of the correlation func- 
tion. 



C{t,T = 0) = Y,Wne 



(16) 



with Wn = |(Ol$|n)|VZ(0). 

However, for any finite temperature (this is always the 
case for finite lattices), all the thermal states with the 
nonzero matrix elements (m|$|n) may contribute to the 
spectral function p{uj). Intuitively in the confinement 
phase, the fundamental degrees of freedom are hadron- 
like modes, thus the thermal states should be multi- 
hadron states. If they interact weakly with each other, 
we can treat them as free particles at the lowest order 
approximation and consider Em as the sum of the ener- 
gies of hadrons including in the thermal state |m). Since 
the contribution of a thermal state |m) to the spectral 
function is weighted by the factor eyip{—Em/T), apart 
from the vacuum state, the maximal value of this fac- 
tor is exp(— Mmin/T) with Mmin the mass of the hght- 
est hadron mode in the system. As far as the quenched 
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glueball system is concerned, the lightest glueball is the 
scalar, whose mass at the low temperature is roughly 
Mg++ 1.6 GeV, which gives a very tiny weight fac- 
tor exp(— Mg++/Tc) 0.003 at Tc in comparison with 
unity factor of the vacuum state. That is to say, for the 
quenched glueballs, up to the critical temperature Tc, 
the contribution of higher spectral components beyond 
the vacuum to the spectral function are much smaller 
than the statistical errors (the relative statistical errors 
of the thermal glueball correlators are always a few per- 
cent) and can be neglected. As a result, the function 
form of p{i-u) in Eq. [15] can be a good approximation for 
the spectral function of glueballs at least up to T^. Ac- 
cordingly, considering the finite extension of the lattice in 
the temporal direction, the function form of the thermal 
correlators can be approximated as 



cosh(M„(l/(2r)-0) 
sinh(Af„/(2r)) ■ 



(17) 



which is surely the commonly used function form for the 
study of hadron masses at low temperatures on the lat- 
tice. As is always done, the glueball masses M„ derived 
by this function are called the pole masses in this work. 



A. Results of the single-cosh fit 

Even though the above discussion are based on the 
weak-interaction approximation for the hadronlike modes 
below Tc, we would like to apply Eq. [17] to analyzing 
the thermal correlators all over the temperature in con- 
cern. The interest of doing so is twofold. First, the ther- 
mal scattering of the glueball-like modes would result in 
a mass shift, say the deviation of the pole mass from 
the glueball mass at zero-temperature, which reflects the 
strength of the interaction at different temperature. Sec- 
ondly, the breakdown of this function form would signal 
the dominance of new degrees of freedom instead of the 
hadronlike modes in the thermal states. 

In practice, after the thermal correlators C{t, T) of the 
optimal operators are obtained according to the steps 
described in Sec. 11(B), the pole masses of the ground 
state (or the lowest spectral component) can be extracted 
straightforwardly. First, for each i?^'-^ channel and at 
each temperature T, the effective mass Mcs{t) as a func- 
tion of t is derived by solving the equation 



C{t + 1, T) _ cosh((t + 1 - Ntl2)atM,a{t)) 



C{t,T) 



cosh((i - Nt/2)atMcs{t)) 



(18) 



Secondly, the effective masses are plotted versus t and 
the plateaus give the fit windows [^1,^2]- Finally, the 
pole masses of the ground states are obtained by fitting 
C{t, T) through a single-cosh function form. As a con- 
vention in this work, we use Mq to represent the mass of 
a glueball state in the physical units and M to represent 
the dimensionless mass parameter in the data processing 
with the relation AI = Meat- 
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FIG. 4: Effective masses at different temperatures in A^^ 
channel. Data points are the effective masses with jackknife 
error bars. The vertical lines indicate the time window [ti, t2] 
over which the single-cosh fittings are carried out, while the 
horizontal lines illustrate the best-fit result of pole masses (in 
each panel the double horizontal lines represent the error band 
estimated by jackknife analysis) 



In Fig.[31 Fig.[Sl Fig.[Sl and Fig.[7]are shown the effec- 
tive masses with jackknife errors at various temperatures 
in Af^, ^r^j and channels, respectively. The 

vertical lines indicate the time window [ti, ^2] over which 
the single-cosh fittings are carried out, while the hori- 
zontal lines illustrate the best-fit result of pole masses 
(in each figure panel the double horizontal lines give the 
error band estimated by jackknife analysis). These fig- 
ures exhibit some common features: At the tempera- 
tures below Tc {Nt = 128, 80, 40), the effective mass 
plateaus show up almost from right the beginning of t, 
as it should be for the optimal glueball operators, while 
at T > Tc {Nt = 36, 24, and 20), the plateaus appear 
later and later in time, and even do not exist at Nt = 20 
(T — 1.90Tc). This observation can be interpreted as 
follows. Since the effective masses are calculated based 
on Eq. I17[ the very early appearance of the plateaus be- 
low Tc implies that the thermal correlators C{t, T) of the 
optimal operators are surely dominated by the ground 
state and can be well described by the function form of 
Eq. [T7] In other words, the picture of weakly interacting 
glueball-like modes makes sense for the state of matter 
below Tc. While at T > Tc, the later appearance and 
the narrower size plateaus signal that the picture of the 
state of matter is distinct from that at T < Tc. How- 
ever, because of the existence of effective mass plateaus 
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FIG. 5: Similar to Fig. [H but in + channel. 



FIG. 7: Similar to Fig.gl but in T^'^ channel. 
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FIG. 6: Similar to Fig. H but in E++ channel. 



up to r - 1.58rc(iVt = 24), the possibility that glueball- 
like modes survive at this high temperature cannot be 
excluded. 

The pole masses in all 20 R^'-^ channels are extracted 
in units a^^ at all temperatures and are shown in Ta- 
ble HVl Specifically, with the lattice spacing determined 
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FIG. 8: The T-dependence of pole masses A~+ 
and glueballs. 
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in Sec. II, the pole masses of A^^ , £■++ and T. 

at T ~ in physical units are M^++=1.576(22)GeV, 
Af^-+=2.488(34)GeV, Me^^ ^ Mj,++=2.364(ll)GeV, 
respectively, which are in agreement with that of previ- 
ous studies |2l| - [28| . From the table, one can see that the 
behaviors of the pole masses with respect to the tempera- 
ture in all 20 channels are uniform: the pole masses keep 
almost constant with the temperature increasing from 
0.30Tc to right below (0.95Tc), and start to reduce 
gradually when T > Tc- When T increases up to l.OOTc, 
the pole masses cannot be extract reliably through the 
single-cosh fit for the lack of clear effective mass plateaus. 
Figure [5] iUustrates these behavior of pole masses in A^^, 
E++ and T++ channels. 
These results imply that glueballs can be very stable 
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TABLE IV: The pole 


xicLSses ( 


in units of ^) 


in all the 20 R^^ 


channels are extracted at all the temperatures. 


dPC 

R 


128 




80 




60 




48 


44 


40 




36 


32 


28 


24 




0.140( 


2) 


0.144( 


3) 


0.144( 


2) 


0.143( 3) 


0.140(2) 


0.140( 


3) 


0.132( 4) 


0.126( 2) 


0.122( 4) 


0.116( 3) 


1 _| 

At 


0.441( 


3) 


0.435( 


3) 


0.434( 


5) 


0.437( 4) 


0.432( 4) 


0.435( 


5) 


0.399( 6) 


0.322( 9) 


0.267(16) 


0.241(13) 




0.221( 


3) 


0.225( 


2) 


0.222( 


2) 


0.225( 2) 


0.218( 3) 


0.222( 


2) 


0.183( 5) 


0.174( 3) 


0.155( 4) 


0.146( 4) 




0.475( 


6) 


0.453( 


8) 


0.447( 


9) 


0.464( 7) 


0.473( 6) 


0.468( 


6) 


0.426(12) 


0.417(10) 


0.287(19) 


0.253(18) 


At+ 


0.323( 


4) 


0.327( 


4) 


0.326( 


4) 


0.330( 2) 


0.326( 4) 


0.332( 


3) 


0.282( 7) 


0.249( 8) 


0.224( 9) 


0.208( 9) 


.J 

At 


0.302( 


5) 


0.308( 


3) 


0.308( 


5) 


0.312( 3) 


0.312( 5) 


0.308( 


6) 


0.268( 6) 


0.241( 7) 


0.220( 8) 


0.201( 6) 


Ao + 


0.450( 


5) 


0.449( 


7) 


0.446( 


5) 


0.440( 6) 


0.452( 4) 


0.448( 


5) 


0.396(10) 


0.340(11) 


0.330(12) 


0.250(14) 


Ar 


0.387( 


3) 


0.388( 


3) 


0.385( 


4) 


0.390( 5) 


0.376( 4) 


0.375( 


4) 


0.354( 7) 


0.293( 7) 


0.268(10) 


0.214( 9) 


£;++ 


0.210( 


1) 


0.205( 


1) 


0.207( 


1) 


0.209(2) 


0.206(1) 


0.189( 


4) 


0.167( 4) 


0.153( 3) 


0.143( 3) 


0.139( 2) 




0.401( 


2) 


0.403( 


2) 


0.401( 


2) 


0.394( 4) 


0.400( 2) 


0.395( 


3) 


0.375( 4) 


0.311( 6) 


0.261( 7) 


0.230( 7) 


E-+ 


0.273( 


1) 


0.266( 


1) 


0.264( 


2) 


0.273( 2) 


0.275( 1) 


0.262( 


2) 


0.218( 4) 


0.196( 4) 


0.183( 4) 


0.181( 4) 


E-- 


0.374( 


1) 


0.368( 


2) 


0.360( 


2) 


0.361( 3) 


0.363( 3) 


0.352( 


4) 


0.308( 8) 


0.262( 6) 


0.231( 6) 


0.213( 6) 


T++ 


0.327( 


2) 


0.326( 


4) 


0.327( 


2) 


0.334( 2) 


0.331( 2) 


0.312( 


5) 


0.287( 7) 


0.266( 3) 


0.227( 6) 


0.215( 4) 


r+- 


0.278( 


1) 


0.274( 


2) 


0.265( 


3) 


0.278( 2) 


0.281( 1) 


0.261( 


3) 


0.207( 6) 


0.199( 2) 


0.181( 4) 


0.175( 2) 




0.372( 


2) 


0.377( 


4) 


0.371( 


3) 


0.380( 2) 


0.374( 2) 


0.370( 


3) 


0.331( 5) 


0.289( 7) 


0.248( 7) 


0.230( 5) 




0.350( 


4) 


0.349( 


2) 


0.344( 


3) 


0.351( 2) 


0.350( 2) 


0.343( 


3) 


0.272( 8) 


0.252( 5) 


0.212( 6) 


0.201( 5) 


rp + + 


0.205( 


1) 


0.209( 


1) 


0.206( 


1) 


0.205( 1) 


0.207( 2) 


0.191( 


3) 


0.160( 3) 


0.152( 2) 


0.148( 2) 


0.143( 2) 




0.322( 


2) 


0.317( 


2) 


0.310( 


4) 


0.317( 3) 


0.320( 2) 


0.303( 


5) 


0.276( 5) 


0.250( 3) 


0.201( 4) 


0.190( 4) 




0.265( 


2) 


0.260( 


3) 


0.264( 


2) 


0.273( 3) 


0.272( 2) 


0.264( 


2) 


0.240( 3) 


0.213( 3) 


0.187( 4) 


0.183( 4) 




0.368( 


2) 


0.358( 


3) 


0.364( 


3) 


0.358( 4) 


0.367( 2) 


0.353( 


5) 


0.282(13) 


0.254( 6) 


0.235( 6) 


0.220( 4) 



below Tc and survive up to l.GTc- This coincides with 
the thermal properties of heavy quarkonia observed by 
model calculation and lattice numerical studies [^-[S lla - 
[TqI . liol l4l| , but different from the observation of a pre- 
vious lattice study on glueballs where the observed pole- 
mass reduction start even at T ~ O.ST^ ISal. 



B. Breit-Wigner analysis 

In the single-cosh analysis, it is seen that, when the 
temperature increases up to Tc, the thermal correlators 
can be well described by Eq. [17] and the pole masses of 
glueballs are insensitive to T. This is in agreement with 
the picture that the state of matter below Tc are made up 
of weakly interacting glueball-like modes. When T > 
the thermal correlators deviate from Eq. [T7] more and 
more. This observation implies that the degrees of free- 
dom are very different from that when T < Tc- Theo- 
retically in the deconfined phase, gluons can be liberated 
from hadrons. However, the study of the equation of 
state shows that the state of the matter right above Tc 
is far from a perturbative gluon gas. In other words, the 
gluons in the intermediate temperature above Tc may 
interact strongly with each other and glueball-like res- 
onances can possibly be formed. Thus different from 
bound states at low temperature, thermal glueballs can 
acquire thermal width due to the thermal scattering be- 
tween strongly interacting gluons and the magnitudes of 
the thermal widths can signal the strength of these types 



of interactions at different temperatures. 

In order to take the thermal width into consideration, 
we also adopt the Breit-Wigner Ansatz, which is sug- 
gested by the pioneering work Ref. [s^, to analyze the 
thermal correlators once more. First, we treat thermal 
glueballs as resonance objects which correspond to the 
poles (denoted by = luq — iT) of the retarded and 
advanced Green functions in the complex oj.plane (note 
that conventionally in particle physics, a resonance pole 
is always denoted as M — iT /2 where M is the mass of 
the resonance and F is its width.) cjq is called the mass 
of the resonance glueball and F its thermal width in this 
work. Secondly, we assume that the spectral function 
p{uj) is dominated by these resonance glueballs. Thus 
the spectral function is parametrized as 



p{uj) = A((5r(w - Wo) — 5t{oj + uoq) + 
where 5, is the Lorentzian function 



5t{x) = — Im 

TT 



1 



1 



(19) 



(20) 



and " . . ." represents the terms of excited states. With 
this spectral function, the thermal glueball correlator 
G{t, T) can be expressed as 



C{t,T) 



duj cosh{uj{^ ~ t)) 
2^ 2sinh(^) 

— OO 

X 27r^ {5r{uj - ujq) - Sr{uj + ojq)) + ... .(21) 
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FIG. 9: Determinations of the fit range [ti,t2] in channel at Nt = 128, 36, and 20. In each row, ui^o°^\t) and F^^^^t) 

obtained by solving Eq. [23] are plotted by data points with jackknife error bars. [ti,t2] are chosen to include the time slices 
between the two vertical lines, where LOg^^^t) and r'"^^(t) show up plateaus simultaneously. The best-fit results of loo and F 
through the function gr{t) are illustrated by the horizontal lines. 



The integral on the right hand side of above equation, denoted by gr{t), can be calculated explicitly as 
I 



grit) = A 



Re 



cosh((a;o + jr)(2^ - t)) 



sinh( 



2T 



oo 

2ujqT ^ cos {2nnTt) 



(27rnT + r)^ 



— (n — > — n) 



(22) 



which can be used as the fit function to extract ujq and F 
from the thermal correlators obtained from the numerical 
calculation. Practically, the infinite series in the above 
equation is truncated by setting the upper limit of the 
summation to be 50, which is tested to be enough for all 
the cases considered in this work. 

In the present study, we carry out the Breit-Wigner 
analysis in A^'^ , ^ -E'^^j a-nd T^~^ channels, whose 
continuum correspondences are ^, and 2++. Al- 
though the variational method is exploited to enhance the 
contribution of the ground state to the thermal correla- 
tors, the contributions from higher spectral components 
cannot be eliminated completely. Therefore, the fit range 
must be chosen properly where the contribution of the 
ground state dominates. We take the strategy advocated 
in Ref.fssj as follows. For a given correlator C{t,T), the 
effective peak position LL)j^^\t) and the effective width 
p(cff)('^^ are obtained by solving the equations 

grjt) ^ C{t,T) 
grit + I) C(t + l,r)' 



gr(^ + l) ^ C{t + l,T) 

grit + 2) Cit + 2,T)' ^ ' 

The statistical errors oiJ^^\t) and F(°ff)(i) can be esti- 
mated through the jackknife analysis. Thus the fit range, 
denoted by [ti , ^2] , is chosen to be the time range where 

u]'^^\t) and T^'^^\t) show up plateaus simultaneously. 
For example, the procedure in T^"*" channel is illustrated 
in Fig. [S] for Nt = 128,36,20 (corresponding to the tem- 
perature T/Tc = 0.30,1.05,1.90), where the fit ranges 
[ti,t2\ are determined to include the time slices between 
the two vertical lines in each figure. 

After the fit ranges for all the thermal correlators are 
chosen, the jackknife analysis can be carried out straight- 
forward and the detailed procedures are omitted here. 
TableEl ED EU and Einishow the fit windows [tiM]. 
the chi-square per degree of freedom ^ jd.o.j , and the 
best-fit results of wq and F at various temperature in 
, A'{^ , E++, and T^^ channels. In almost aU the 
cases, the fit ranges start from ti = 1, 2, or 3, and last for 
quite a few time slices. This reflects that, as is expected, 
the optimal glueball operators couple almost exclusively 
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TABLE V: The best-fit ujq and F of A++ channel at different 
T through the Breit-Wigner fit. Also listed are the fit window 
[ti,f2] and the chi-square per degree of freedom, jd.o.f . 



Nt 


T/Tc 




F 




x'/d.o.f 


128 


0.30 


0.142(2) 


0.008(10) 


(3, 6) 


1.266 


80 


0.47 


0.146(2) 


0.013(9) 


(2, 4) 


0.921 


60 


0.63 


0.144(2) 


0.008(3) 


(2, 7) 


0.135 


48 


0.79 


0.143(2) 


0.014(6) 


(2, 4) 


0.639 


44 


0.86 


0.142(2) 


0.004(3) 


(1, 7) 


0.758 


40 


0.95 


0.143(2) 


0.003(4) 


(4, 7) 


0.850 


36 


1.05 


0.146(2) 


0.028(4) 


(3, 6) 


0.960 


32 


1.19 


0.141(2) 


0.053(2) 


(1,4) 


0.393 


28 


1.36 


0.142(2) 


0.056(4) 


(2, 5) 


1.253 


24 


1.58 


0.143(1) 


0.059(3) 


(1, 4) 


0.302 


20 


1.90 


0.146(2) 


0.077(5) 


(2, 4) 


0.918 
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FIG. 11: F's are plotted versus T/Tc for A-+ E++ , 

and T^^ channels. The vertical lines indicate the critical 
temperature. 



TABLE VI: The best-fit a;o and T of A~+ channel at different 
T through the Breit-Wigner fit. Also listed are the fit window 
[ti,t2] and the chi-square per degree of freedom, /d.o.f. 



Nt 


T/T, 




F 




xVd.o.f 


128 


0.30 


0.226(2) 


0.006(5) 


(3, 9) 


0.509 


80 


0.47 


0.228(2) 


0.008(4) 


(2, 6) 


0.640 


60 


0.63 


0.227(1) 


0.013(5) 


(2, 7) 


0.216 


48 


0.79 


0.228(2) 


0.012(5) 


(2, 6) 


0.177 


44 


0.86 


0.229(2) 


0.013(4) 


(2, 8) 


0.184 


40 


0.95 


0.224(2) 


0.004(6) 


(3, 6) 


0.549 


36 


1.05 


0.221(2) 


0.037(3) 


(1,8) 


0.935 


32 


1.19 


0.219(2) 


0.047(3) 


(1,6) 


0.250 


28 


1.36 


0.211(3) 


0.068(6) 


(2, 5) 


0.091 


24 


1.58 


0.208(3) 


0.089(7) 


(2, 4) 


0.003 


20 


1.90 


0.211(3) 


0.099(9) 


(2, 6) 


0.083 



to the lowest spectral components after the implemen- 
tation of the variational method. All the /d.o.f^s are 
^ 0(1) or even smaller, which reflect the reliability of 
the fits. 

The main features of the best fit wq and F based on 
Breit-Wigner Ansatz are described as follows: 

• The peak positions loq of the spectral functions p{uj) 
are insensitive to the temperature in all the consid- 
ered channels. In particular, the wq in A\^~^ chan- 
nel keeps almost constant all over the temperature 
range from O.SOTc to 1.90Tc. In the other three 
channels, the wo's do not change within errors be- 
low Tc, but reduce mildly with the increasing tem- 
perature above Tc- The reduction of wq at the high- 
est temperature T — 1.90Tc is less than 5% in these 
three channels. 

• In all four channels, the thermal widths F are small 
and do not vary much below Tc, but grow rapidly 
with the increasing temperature when T > Tc- Be- 
low Tc, the thermal widths are of order F ~ 5% or 
even smaller (especially for the A^^ F is consistent 
with zero). The thermal widths increase abruptly 



when the temperature passes Tc and reach values 
- cjo/2 at T = 1.90Tc. 

These features can be seen easily in Fig. [TU] and [TTl 
where the behaviors of luq and F with respect to the 
temperature T are plotted for all four channels. The 
line_shapes of the spectral functions with the best-fit pa- 
rameters at different T are shown in Fig. [12] for 
-E"*""*", and T^^ channels (we do not plot the spec- 
tral function of A^^ channel due to the small thermal 
widths). 



IV. SUMMARY AND DISCUSSIONS 

On 24"^ X Nt anisotropic lattices with the anisotropy 
^ = 5 at the gauge coupling /3 = 3.2, the thermal 
glueball correlators are calculated in a large tempera- 
ture range from O.SOTc to 1.90Tc, which are realized by 
varying Nt to represent different temperatures. Based 
on the lattice spacing as = 0.0878(4) fm determined by 
Tg'^ — (410(20) MeV), the spatial extension of the lattices 



11 



50 
45 
40 
35 
30 
25 
20 
15 
10 
5 




(a) A- 





T/T =1.58 




T/Tj=1.05 




T/T =0.95 




T/T =0.30 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 




(c) T++ 



T/T =1.58 
T/Tj=1.05 
T/Tj=0.95 
T/T =0.30 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



FIG. 12: Plotted are the spectral function p{ui) at T /Tc — 0.30, 0.95,1.05 and 1.58 with the best-fit parameters. Panel (a)_(c) 
are for A^^ , and T^^ channels, respectively. 



TABLE VII: The best-fit ljq and F of E++ channel at different 
T through the Breit-Wigner fit. Also listed are the fit window 
[ii,i2] and the chi-squarc per degree of freedom, /d.o.f. 



Nt 




[Jo 


F 




xVdof 


128 


0.30 


0.212(1) 


0.012(4) 


(2, 5) 


0.274 


80 


0.47 


0.211(1) 


0.006(3) 


(2, 8) 


0.616 


60 


0.63 


0.212(1) 


0.010(3) 


(2, 9) 


0.844 


48 


0.79 


0.213(1) 


0.011(3) 


(2, 5) 


0.206 


44 


0.86 


0.211(1) 


0.011(3) 


(2, 8) 


1.268 


40 


0.95 


0.207(1) 


0.022(3) 


(2, 6) 


0.250 


36 


1.05 


0.205(2) 


0.034(2) 


(1,8) 


0.183 


32 


1.19 


0.200(1) 


0.049(2) 


(1,6) 


0.478 


28 


1.36 


0.191(2) 


0.067(4) 


(2, 6) 


0.297 


24 


1.58 


0.189(2) 


0.083(5) 


(2, 4) 


0.253 


20 


1.90 


0.196(2) 


0.091(5) 


(2, 4) 


0.046 



TABLE VIII: The best-fit wo and F of T++ channel at dif- 
ferent T through the Breit-Wigner fit. Also listed are the fit 
window [ti,t2] and x^/d-o.f. 



Nt 






F 






128 


0.30 


0.210(1) 


0.008(2) 


(3,10) 


0.442 


80 


0.47 


0.213(1) 


0.009(3) 


(3, 9) 


0.696 


60 


0.63 


0.213(1) 


0.012(3) 


(3, 7) 


0.326 


48 


0.79 


0.210(1) 


0.007(4) 


(3, 6) 


0.288 


44 


0.86 


0.214(1) 


0.012(2) 


(2, 6) 


0.437 


40 


0.95 


0.208(1) 


0.023(1) 


(1,7) 


0.743 


36 


1.05 


0.204(1) 


0.039(2) 


(1,7) 


0.606 


32 


1.19 


0.199(1) 


0.047(1) 


(1,6) 


0.527 


28 


1.36 


0.196(2) 


0.064(3) 


(2, 5) 


0.022 


24 


1.58 


0.194(1) 


0.077(4) 


(2, 7) 


0.119 


20 


1.90 


0.196(2) 


0.093(4) 


(2, 5) 


0.027 



are estimated to be (2.1 fm)'^, which is large enough to 
be free of the finite volume effects. On the other hand, 
because of the large anisotropy, there are enough data 
point in the temporal direction for the thermal correla- 
tors to be analyzed comfortably even at the highest tem- 
perature T ~ 2Tc concerned in this work. With the im- 
plementation of the smearing scheme and the variational 



method, we can construct the optimal glueball operators 
in all the symmetry channel, which couple mostly to the 
lowest-lying states (or more precisely, the lowest-lying 
spectral components). As a result, the thermal corre- 
lators of these operators can be considered to be con- 
tributed dominantly from these lowest-lying states. The 
thermal correlators are analyzed based on two ansatz, 
say, the single-cosh function form and the Breit-Wigner 
ansatz. In Table |TX] and Table H the "pole masses 
Mg obtained by single-cosh analysis and (wq, F) obtained 
based on the Breit-Wigner ansatz are combined together 
for comparison (all the data are converted into physical 
units). 

The most striking observation from the single-cosh 
analysis is that, in all 20 R^'^ channels, the best-fit pole- 
masses Mg are almost constant within errors from the 
low temperature up to right below the critical tempera- 
ture Tc. This is what should be from the point of view of 
deconfinement phase transition of QCD: Since below 
the system is in the confinement phase, the fundamen- 
tal degrees of freedom must be hadrons. Above Tc, the 
reduction of the pole masses does signal the QCD transi- 
tion, after which the state of the matter is very different 
from that below Tc. However, the existence of effective 
mass plateaus, from which the pole masses are extracted, 
also implies that color singlet objects, the glueball-like 
modes, can also survive at the intermediate temperature 
above Tc. The results of the Breit-Wigner fit are consis- 
tent with this picture. In the Breit-Wigner ansatz, ther- 
mal widths F are introduced to glueball states to account 
for the effects of finite temperature, such as the thermal 
scattering and the thermal fluctuations. As shown in Ta- 
ble [IX] and |Xl below Tc (or in the confinement phase), 
the best-fit wq's are very close to the pole masses, and 
the thermal widths F are very tiny and are always of a 
few percent of wq. This means the glueball states are 
surely stable in the confinement phase and the thermal 
interaction among them are weak. With the temperature 
increasing above Tc, while the temperature dependence of 
wq's is very mild, the thermal widths F grow rapidly and 
reach values of roughly half of w's at T ^ 1.9Tc. This 
clearly reflects that glueballs act as resonances are un- 
stable more and more, and the reduction of pole masses 
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TABLE IX: The "pole masses" Ma obtained by single-cosh analysis and {uJo,V) obtained based on the Breit-Wigner ansatz 
are combined together for comparison. Listed in the table are the results in A^^ and A'l'^ channels (all the data are converted 
into the physical units). 





AT* 


Nt 




mcpeV] 


a;o[GeV] 


r[GeV] 


mG[GeV] 


a;o[GeV] 


r[GeV] 


128 


0.30 


1.576(22) 


1.602(14) 


0.091(113) 


2.488(31) 


2.549(17) 


0.069(52) 


80 


0.47 


1.621(29) 


1.644(22) 


0.145(100) 


2.533(24) 


2.560(17) 


0.086(44) 


60 


0.63 


1.627(18) 


1.621(16) 


0.098(38) 


2.499(27) 


2.559(16) 


0.147(54) 


48 


0.79 


1.616(21) 


1.612(26) 


0.156(67) 


2.533(23) 


2.564(20) 


0.135(55) 


44 


0.86 


1.577(25) 


1.598(17) 


0.045(34) 


2.454(34) 


2.577(18) 


0.144(48) 


40 


0.95 


1.576(39) 


1.621(20) 


0.034(46) 


2.499(25) 


2.525(26) 


0.042(71) 


36 


1.05 


1.486(43) 


1.638(28) 


0.315(48) 


2.060(48) 


2.490(23) 


0.413(32) 


32 


1.19 


1.418(21) 


1.588(25) 


0.586(28) 


1.959(37) 


2.464(20) 


0.529(28) 


28 


1.36 


1.373(48) 


1.599(26) 


0.619(43) 


1.745(43) 


2.375(28) 


0.768(71) 


24 


1.58 


1.306(41) 


1.613(28) 


0.664(37) 


1.644(45) 


2.308(32) 


0.999(83) 


20 


1.90 




1.642(32) 


0.873(59) 




2.380(35) 


1.114(99) 



TABLE X: The pole masses Mq obtained by single-cosh analysis and {u)o,V) obtained based on the Breit-Wigner ansatz are 
combined together for comparison. Listed in the table are the results in E'^^ and T^"*" channels (all the data are converted 
into physical units). 



E++ 


rn+ + 


Nt 


T/T, 


mG[GeV] 


a;o[GeV] 


r[GeV] 


mG[GeV] 


ajo[GeV] 


r[GeV] 


128 


0.30 


2.364(11) 


2.385(12) 


0.140(42) 


2.308(14) 


2.363(10) 


0.091(25) 


80 


0.47 


2.308(13) 


2.368(12) 


0.069(29) 


2.353(15) 


2.387(10) 


0.105(32) 


60 


0.63 


2.330(11) 


2.383(12) 


0.116(32) 


2.319(13) 


2.396(10) 


0.140(32) 


48 


0.79 


2.353(19) 


2.393(15) 


0.129(34) 


2.308(15) 


2.362(14) 


0.083(43) 


44 


0.86 


2.319(15) 


2.379(12) 


0.119(32) 


2.330(21) 


2.405(10) 


0.136(26) 


40 


0.95 


2.263(14) 


2.327(15) 


0.247(38) 


2.218(16) 


2.344(11) 


0.259(16) 


36 


1.05 


1.880(41) 


2.305(17) 


0.382(23) 


1.925(33) 


2.298(14) 


0.437(17) 


32 


1.19 


1.722(35) 


2.247(16) 


0.549(20) 


1.801(25) 


2.244(11) 


0.532(15) 


28 


1.36 


1.610(31) 


2.155(19) 


0.754(43) 


1.666(23) 


2.205(17) 


0.717(36) 


24 


1.58 


1.565(23) 


2.132(21) 


0.937(55) 


1.610(27) 


2.184(16) 


0.870(41) 


20 


1.90 




2.201(23) 


1.023(60) 




2.209(20) 


0.935(48) 



above Tc can be taken as the effect of these growing ther- 
mal widths. 

To summarize, in pure gauge theory, the state of mat- 
ter is dominated by weakly interacting hadronlike states 
below Tc] when T > T^, glueball states survive as res- 
onancelike modes up to a temperature T ~ 1.9Tc with 
their thermal widths growing with increasing T, which 
implies that in this intermediate temperature range, glue- 
balls are unstable and may decay into gluons, and re- 
versely gluons also interact strongly enough to form 
glueball-like resonances. The two procedure may reach 
the thermal equilibrium at a given temperature, such 
that the gluon degree of freedom become more and more 
important with T increasing. At very high temperature, 
the glueball-like resonances may disappear finally and the 
state of matter can thereby be described by a perturba- 
tive gluon plasma. This picture is coincident with the ob- 
servations both in the study of equation of state of QCD 
and the thermal properties of heavy quarkonia. On the 
other hand, the surprising results of RHIC experiments 
mav also suDDort this oicture to some extent. First, the 



data of RHIC experiments are well described by the hy- 
drodynamical model[42,]. Secondly, the investigation of 
elliptic flow data using a Boltzmann-type equation for 
gluon scattering is not consistent with the perturbative 
QCD apparently [i^. So the quark-gluon plasma at the 
RHIC temperature is most likely a strongly interacting 
system. 
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Thermal properties of glueballs in SU (3) Yang- Mills theory are investigated in a large temperature 
range from 0.3Tc to l.QTc on anisotropic lattices. The glueball operators are optimized for the 
projection of the ground states by the variational method with a smearing scheme. Their thermal 
correlators are calculated in all 20 symmetry channels. It is found in all channels that the pole 
masses Mq of glueballs remain almost constant when the temperature is approaching the critical 
temperature from below, and start to reduce gradually with the temperature going above Tc- The 
correlators in the 0^"*", and 2^"*" channels are also analyzed based on the Breit-Wigner Ansatz 
by assuming a thermal width F to the pole mass ojq of each thermal glueball ground state. While 
the values of ujo are insensitive to T in the whole temperature range, the thermal widths F exhibit 
distinct behaviors at temperatures below and above Tc. The widths are very small (approximately 
few percent of ujq or even smaller) when T < Tc, but grow abruptly when T > Tc and reach values 
of roughly F ~ a;o/2 at T 1.9Tc. 

PACS numbers: 12.38.Gc, 11. 15. Ha, 14.40.Rt, 25.75.Nq 



I. INTRODUCTION 

The past two or three decades witnessed intensive and 
extensive studies on the t)hase transition of quantum 
chroniodynamics(QCD) [J], which is believed to be the 
fundamental theory of strong interaction. Based on the 
two characteristics of QCD, namely the conjectured color 
confinement at low energies and the asymptotic freedom 
of gluons and quarks at high energies, QCD at finite tem- 
perature is usually described by two extreme pictures. 
One is with the weakly interacting meson gas in the 
low temperature regime and another is with perturba- 
tive quark gluon plasma (QGP) in the high tempera- 
ture regime. The two regimes are bridged by a decon- 
flnement phase transition (or crossover). The study of 
the equation of state shows that the perturbative pic- 
ture of QGP can only be achieved at very high temper- 
atures T > 2Tc. In other words, the dynamical degrees 
of freedom up to the temperature of a few times of Tc 
are not just the quasifree gluons and quarks 0]. Some 
other theoretical studies also support this scenario and 
conjecture that in the intermediate temperature range 
above there may exist different types of excitations 
corresponding to different distance scales [1, |3| render- 
ing the thermal states much more complicated. Apart 
from the quasifree quarks and gluons at the small dis- 
tance scale, the large scale excitations can be effective 
low-energy modes in the mesonic channels as a result 
of the strongly interacting partons Q . The properties of 
the interaction among quarks and gluons at low and high 
temperatures can be studied with thermal correlators. 



There have been many works on the correlators of char- 
monia at finite temperature. Phenomenological studies 
predicted the binding between quarks is reduced to dis- 
solve J/tp at temperatures close to Tc and proposed the 
suppression of charmonia as a signal of QGP [6, 7]. For 
example, potential model studies show that excited states 
like Tp and Xc are dissociated at Tc, while the ground 
state charmonia J/^ and r]c survive up to T = l.lTc [1- 
[isj . However, it is unclear whether the potential model 
works well at finite temperatures p^ . In contrast, many 
recent numerical studies indicate that J/^* and rjc might 
still survive above 1.5Tc p"5l - [l9j . Of course, it is pos- 
sible that the cc states observed in lattice QCD are 
just scattering states. A further lattice study on spa- 
tial boundary-condition dependence of the energy of low- 
lying cc system concludes that they are spatially local- 
ized (quasi)bound states in the temperature region of 
1.11 ~ 2.07Tc [l^]. Obviously, the results of numerical 
lattice QCD studies are coincident to the picture of the 
QCD transition in the intermediate temperature regime. 

Until now most of the lattice studies on hadronic corre- 
lators are in the quenched approximation. Because of the 
lack of dynamical quarks in quenched QCD the binding 
of quark-antiquark systems must be totally attributed 
to the nonperturbative properties of gluons, which are 
the unique dynamical degree of freedom in the theory. 
Since glueballs are the bound states of gluons, a natural 
question is how glueballs respond to the varying tem- 
peratures. At low temperature T ^ 0, the existence 
of quenched glueballs have been verified by extensive 
lattice numerical studies, and their spectrum are also 
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established quite well |21| - |28l |. An investigation of the 
evolution of glueballs versus the increasing temperature 
is important to understand the QCD transition [2^ [s^l 
and the hadronization of quark-gluon plasma (3l| . From 
the point of view of QCD sum rules, glueball masses are 
closely related to the gluon condensate. Lattice stud- 
ies [12] and model calculations indicate that the 
gluon condensate keeps almost constant below Tc and 
reduces gradually with the increasing temperature above 
Tc- Based on this picture, it is expected intuitively that 
glueball masses should show a similar behavior also until 
they melt into gluons [S^l- In fact, there has already been 
a lattice study on the scalar and tensor glueball proper- 
ties at finite temperature [ssj . In contrast to the expecta- 
tion and the finite T behavior of charmonium spectrum, 
it is interestingly observed that the pole-mass reduction 
starts even below Tc (mG(T ~ Tc) ~ 0.8tog(T - 0)). It 
is known that the spatial symmetry group on the lattice 
is the 24-element cubic point group O, whose irreducible 
representations are R = Ai, A2, E, Ti, and T2. Along 
with the parity P and charge conjugate transformation 
C, all the possible quantum numbers that glueballs can 
catch are R^'^ with PC = ++, — h, H — ,, and -I— 1-, which 
add up to 20 symmetry channels. Motivated by the dif- 
ferent temperature behaviors of cc systems with different 
quantum numbers, we would like to investigate the tem- 
perature dependence of glueballs in this paper. 

Our numerical study in this work is carried out on 
anisotropic lattices with much finer lattice in the tem- 
poral direction than in spatial ones. In order to explore 
the temperature evolution of glueball spectrum, the tem- 
perature range studied here extends from 0.3Tc to 1.9Tc, 
which is realized by varying the temporal extension of 
the lattice. Using anisotropic lattices, the lattice param- 
eters are carefully determined so that there are enough 
time slices for a reliable data analysis even at the high- 
est temperature. In the present study, we are only inter- 
ested in the ground state in each symmetry channel R^'~^ . 
For the study optimized glueball operators that couple 
mostly to the ground states are desired. Practically, these 
optimized operators are built up by the combination of 
smearing schemes and the variational method [2l| - |23| . 
In the data processing, the correlators of these optimized 
operators are analyzed through two approaches. First, 
the thermal masses Mq of glueballs are extracted in all 
the channels and all over the temperature range by fit- 
ting the correlators with a single-cosh function form, as is 
done in the standard hadron mass measurements. Thus 
the T-cvolution of the thermal glueball spectrums is ob- 
tained. Secondly, with respect that the finite tempera- 
ture effects may result in mass shifts and thermal widths 
of glueballs, we also analyze the correlators in A'^'^ , ^ 
E++, and T^'^ channels with the Breit-Wi gner Ansatz 
which assumes these glueball thermal widths, say, change 
Mq into ljJq — iT in the spectral function (see below). As 
a result, the temperature dependence of wq and T can 
shed some light on the scenario of the QCD transition. 

This paper is organized as follows. In Sec. II, a descrip- 



tion of the determination of working parameters, such as 
the critical temperature Tc, temperature range, and lat- 
tice spacing a^, as well as a brief introduction to the vari- 
ational method is given. In Sec. Ill, after a discussion of 
its feasibility, the results of the single-cosh fit to the ther- 
mal correlators are described in details. The procedure of 
the Breit-Wigner fit is also given in this section. Section 
IV gives the conclusion and some further discussions. 



II. NUMERICAL DETAILS 

For heavy particles such as charmonia and glueballs, 
the implementation of anisotropic lattices is found to be 
very efficient in the previous numerical lattice QCD stud- 
ies both at low and finite temperatures. On the other 
hand, the Symanzik improvement and tadpole improve- 
ment schemes of the gauge action are verified to have bet- 
ter continuum extrapolation behaviors for many physical 
quantities. In other words, the finite lattice spacing ar- 
tifacts are substantially reduced by these improvements. 
With these facts, we adopt the following improved gauge 
action which has been extensively used in the study of 
glueballs [Hl-il, 
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where /? is related to the bare QCD coupling constant, 
^ = as/ at is the aspect ratio for anisotropy (we take ^ = 5 
in this work), Us and ut are the tadpole improvement 
parameters of spatial and temporal gauge links, respec- 
tively. rJc = Ec \Re-Tr{\ - Wc), with Wc denoting 
the path-ordered product of link variables along a closed 
contour C on the lattice, fi^p includes the sum over all 
spatial plaquettes on the lattice, ^Itp includes the tem- 
poral plaquettes , VLgr denotes the product of link vari- 
ables about planar 2x1 spatial rectangular loops, and 
fistr refers to the short temporal rectangles (one tempo- 
ral link, two spatial). Practically, Ut is set to 1, and Ug is 
defined by the expectation value of the spatial plaquette, 
\TrP^^, >i/4. 



A. Determination of critical temperature 

Since the temperature T on the lattice is defined by 

1 



Ntat' 



(2) 



where Nt is the temporal lattice size, T can be changed 
by varying either Nt or the coupling constant /3 which 
is related directly to the lattice spacing. In order for 
the critical temperature to be determined with enough 
precision, for a given Nt — 24, we first determine the 
critical coupling 13c, because j3 can be changed continu- 
ously. The order parameter is chosen as the susceptibility 
Xp of Polyakov line, which is defined as 



(3) 
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TABLE I: The simulation parameters for the determination 
of the critical point. The configurations are selected every ten 
sweeps. 

f} Total configurations Thermalization Bin size 



2.80 
2.805 
2.81 
2.815 
2.82 



20000 
30000 
30000 
20000 
8000 



5000 
10000 
10000 
5000 
3000 



1000 
1000 
1000 
1000 
500 




2.S00 2,802 2,804 2,806 2,808 2,810 2,812 2,814 2,816 2,8 



FIG. 1: The xp extrapolation based on the spectral density 
method. The open triangles denote the simulated values of 
Xp, while the filled squares are the extrapolated values. The 
peak position gives the critical /3c ~ 2.808. 



where Q is the Z{3) rotated Polyakov line, 

RePexp[-27ri/3]; arg P e [tt/B, tt) 
e = i ReP; argPe [-7r/3,7r/3) ,(4) 

RePexp[27rV3]; arg P e [-tt, -7r/3) 

and P represents the trace of the spatially averaged 
Polyakov line on each gauge configuration. 

After a /3-scanning on = 24"* anisotropic lattices 
with ^ 5, the critical point is trapped in a very nar- 
row window G [2.800, 2.820]. In order to determine T^, 
more precisely, a more refined study is carried out in the 
/3 window mentioned above with much larger statistics 
through the spectral density method. Practically, the 
spectral density method [3B,|37| is applied to extrapolate 
the simulated xp's at /3 = 2.805, 2.810, and 2.815. In ta- 
bleUare the numbers of heat-bath sweeps for each f3. The 
extrapolation results are illustrated in Fig. [1] where the 
open triangles denote the simulated values of xpi while 
the filled squares are the extrapolated values. Finally, 
the peak position gives the critical coupling constant 
(3c = 2.808, which corresponds to the critical tempera- 
ture Tc « 0.724rn ^ = 296 MeV with the lattice spacing 
ro/a, = 3.476 HI and r^^ = 410(20) MeV. 

With Tc fairly determined, the working coupling con- 
stant /? is set based on two requirements. First, the spa- 
tial volume of the lattice should be large enough in order 
for the glueballs to be free of any sizable finite volume 
effects. Secondly, we require that temporal lattice has a 



TABLE II: Listed are the parameters used to check the critical 
behavior for /3=3.2. The configurations are selected every ten 
sweeps. 

Nt Total configurations Thermalization < P > xp 



60 
48 
44 
40 
36 
32 
28 
24 
20 



2000 
2000 
8000 
8000 
8000 
3000 
2000 
2000 
2000 



500 


-8.73x10- 


-5 


6.65x10 


500 


6.01x10" 


5 


1.81x10 


2000 


2.25x10" 


3 


3.12x10 


2000 


1.72x10" 


2 


9.14x10 


2000 


5.21x10" 


2 


3.10x10 


1000 


8.51x10" 


2 


2.23x10 


500 


0.1253 




2.00x10 


500 


0.1817 




2.09x10 


500 


0.2571 




1.82x10 



0.014 
0.012 
0.01 
0.008 
0.006 
0.004 
0.002 
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20 25 30 35 40 45 50 55 60 
N,=l/(Ta,) 

FIG. 2; Xp is plotted versus Nt a.t (3 — 3.2. There is a peak 
of Xp near Nt — 40. 



good resolution even at the temperature T ^ 2Tc. Prac- 
tically the working coupling constant is finally set to be 
/3 = 3.2. The lattice spacing at this (3 is set by calcu- 
lating the static potential V{r) on an anisotropic lattice 
24^^ x 128. With the conventional parametrization of V(r), 



V{r) = Vo + ar + 



(5) 



the lattice spacing Og is determined in the units of ro to 
be 



1.65- 



0.1825(7) 



(6) 



where ro is the hadronic scale parameter. If we take 
= 410(20) Afey, we have = 0.0878(4) fm. The 
spatial volume at L = 24 is therefore estimated to be 
(2.1fm)3. On the other hand, using = 296MeV ob- 
tained at (3 = 2.808 as a rough estimate of and ignor- 
ing the systematic error due to finite lattice spacings, Tc 
and 2Tc at (3 = 3.2 are expected to be achieved around 
Nt ^ 40 and Nt ~ 20, respectively. Obviously, the above 
two requirements are all satisfied. 

Based on the discussions above, with a fixed (3 = 3.2, 
the calculations of the thermal correlators of glueballs are 
carried out on a series of lattice 24^ x Nt with Nt = 20, 



4 



r 




FIG. 3: Prototype Wilson loops used in making the smeared 
glueball operators [2^. 



24, 28, 32, 36, 40, 44, 48, 60, 80, and 128, which cover the 
temperature range 0.3Tc < T < 2Tc. As a cross-check, 
Xp at different Nt are calculated first and the results are 
shown in Fig. [2] and Table |lTl It is clear that the ex- 
pectation value of the Polyakov line drops to zero near 
Nt = 40 and the peak position of xP: which gives the 
critical temperature, is trapped between Nt — 36 and 
Nt — AO. In practice, we do not carry out a precise de- 
termination of Tc at /3 = 3.2, but take the temperature 
at Nt = 38, T w (38at)-i = (38a^/0"^ = 296MeV, as 
an approximation of Tc (/3 — 3.2), to scale the temper- 
atures involved in this work. It should be noted that, 
owing to the lattice artifact, the critical temperature Tc 



determined at different lattice spacing (or (3) may differ 
from each other. The closeness of Tc(/3 = 2.808) and 
Tc{P = 3.2) may signal that the lattice spacing depen- 
dence of Tc is mild in this work due to the application of 
the improved gauge action. 



B. Variational method 

It is known that many states contribute to a hadronic 
two-point function. Ideally one can extract the informa- 
tion of the lowest-lying states from the two-point function 
in the large time region if it lasts long enough in the time 
direction. This is the case for some light hadron states, 
such as TT meson, K meson, etc. However, for heavy par- 
ticles, especially for glueballs whose correlation function 
are much more noisy than that of conventional hadrons 
made up of quarks, their two-point functions damp so 
fast with time that they are always undermined by noise 
rapidly before the ground states dominate. Practically, 
in the study of the glueball sector, in order to enhance the 
overlap of the glueball operators to the ground state, the 
commonly used techniques are the smearing schemes and 
the variational techniques. In this work, we adopt the 
sophisticated strategy implemented by the studies of the 
zero-temperature glueball spectrum [2l| - |23j , which is out- 
lined below. First, for each gauge configuration, we per- 
form six smearing/fuzzing schemes to the spatial links, 
which are various combinations of the single-link proce- 
dure (smearing) and the double-link procedure (fuzzing) 



U-^ix) = Psu(3){Ujix) + Xs Uk{x)U,ix + k)ul{x + j)}, 

±(Mj) 

Ujix) = Psu{3){UAx)U,{x+j) + Xf Uk{x)U, {x + k)U,{x+ J + k)Uk{x + 2j)}, (7) 

I 



where Psu(3) denotes the projection into SU{3) and is 
realized by the Jacobi method [s^. The six schemes are 
given explicitly as s\°, s\^^, sf^, hf<S)si°^ hf<S)s\l, 
f\f ^ where s/ f denotes the smearing/fuzzing pro- 
cedure defined in Eq. (O, and Xs/Xf the tunable pa- 
rameter which we take As = 0.1 and A/ — 0.5 in this 
work. Secondly, we choose the same prototype Wilson 
loops as that in Ref. [1^ (as shown in Fig. such 
that for each smearing/fuzzing scheme, all the different 
spatially oriented copies of these prototypes are calcu- 
lated from the smeared gauge configurations. Thus for 
a given irreducible representation R of the spatial sym- 
metry group O, say, R = Ai, A2, E,Ti, or T2, a real- 
ization of R can be a specific combination of differently 
oriented Wilson loops generated from the same proto- 
type loop (one can refer to Ref. [23l| for the concrete 



combinational coefficients). The glueball operators (j) 
with the quantum number R^'~^ are thereby constructed 
along with the spatial refiection and the time inversion 
operations. In practice, we establish four realizations 
of each R^'~^ which are based on four different proto- 
types, respectively. Therefore, along with the six smear- 
ing/fuzzing schemes, an operator set of the same specific 
quantum number i?^*^ is composed of 24 different op- 
erators, {0Q, a — 1,2,..., 24}. The last step is the im- 
plementation of the variational method (VM). The main 
goal of VM is to find an optimal combination of the set of 
operators, •l' = ^a'/'ai which overlaps most to a specific 
state (in this work, we only focus on the ground states). 
The combinational coefficients v = {va, a = 1, 2, . . . , n} 



5 



can be obtained by minimizing the effective mass, 

m{tD) = -—\n^ . , (8) 

at ti) — I, where C'afj (t) is the correlation matrix of the 
operator set, 



^a^W = I](O|0a(i + r)<^^(r)|O). 



(9) 



This is equivalent to solving the generalized eigenvalue 
equation 



(10) 



and the eigenvector v gives the desired combinational co- 
efficients. Thus, the optimal operator that couples most 
to a specific states (the ground state in this work) can be 
built up as 



$ = ^ Va(j)a, 



(11) 



whose correlator C{t) is expected to be dominated by the 
contribution of this state. 



III. DATA ANALYSIS OF THE THERMAL 
CORRELATORS OF GLUEBALLS 

AU 20 R^^' channels, with R = Ai, A2, E,Ti,T2 

and PC = ++,H — +, , are considered in the 

calculation of the thermal correlators of glueballs on 
anisotropic lattices mentioned in Sec. II. At each tem- 
perature, after 10000 pseudo-heat-bath sweeps of ther- 
malization, the measurements are carried out every three 
compound sweeps, with each compound sweep composed 
of one pseudo-heat-bath and five micro-canonical over- 
relaxation(OR) sweeps. In order to reduce the possible 
autocorrelations, the measured data are divided into bins 
of the size rimb = 400, and each bin is regarded as an in- 
dependent measurement in the data analysis procedure. 
The numbers of bins A^bin and n,„b at various tempera- 
ture are listed in Table Hill 

Theoretically, under the periodic boundary condi- 
tion in the temporal direction, the temporal correlators 
C{t, T) at the temperature T can be written in the spec- 
tral representation as 



^ |(n|$|m)|2 f En 
= ^TTTTTTT^f^ exp 



2Z{T) 



En 



2T 



X cosh 



dujp{uj)K(uj, T), 



(12) 



TART P, TTT 


: Simulation parameters to calculate 


ilLlcUaii nJJCU 


■fniTtl r? — 

LI Llill. LJ 


3.2, as = 0.0878 fm, Ls 


= 2.11 fm. 








?^mb 




128 


0.30 


400 


24 


80 


0.47 


400 


30 


60 


0.63 


400 


44 


48 


0.79 


400 


40 


44 


0.86 


400 


44 


40 


0.95 


400 


40 


36 


1.05 


400 


40 


32 


1.19 


400 


56 


28 


1.36 


400 


40 


24 


1.58 


400 


40 


20 


1.90 


400 


40 



with a T-dcpendent kernel 

cosh(w/(2T) - wt) 



sinh(w/(2r)) 



and the spectral function, 
^(-) = L 2Z{T) ' 

X (,5(W - {En - Era) - Siuj - {En, - S„)),(14) 

where Z{T) is the partition function at T, and £"„ the 
energy of the thermal state \n) (|0) represents the vac- 
uum state). In the zero-temperature limit (T — 0), due 
to the factor exp(— £^„i/r), the spectral function p{lli) 
degenerates to 



p(^) = E' 



2Z{0) 



{S{uj - En) - S{u; + En)) , (15) 



thus we have the function form of the correlation func- 
tion. 



C{t,T = 0) = Y,Wne 



(16) 



with Wn = |(Ol$|n)|VZ(0). 

However, for any finite temperature (this is always the 
case for finite lattices), all the thermal states with the 
nonzero matrix elements (m|$|n) may contribute to the 
spectral function p{uj). Intuitively in the confinement 
phase, the fundamental degrees of freedom are hadron- 
like modes, thus the thermal states should be multi- 
hadron states. If they interact weakly with each other, 
we can treat them as free particles at the lowest order 
approximation and consider Em as the sum of the ener- 
gies of hadrons including in the thermal state |m). Since 
the contribution of a thermal state |m) to the spectral 
function is weighted by the factor eyip{—Em/T), apart 
from the vacuum state, the maximal value of this fac- 
tor is exp(— Mmin/T) with Mmin the mass of the hght- 
est hadron mode in the system. As far as the quenched 
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glueball system is concerned, the lightest glueball is the 
scalar, whose mass at the low temperature is roughly 
Mg++ 1.6 GeV, which gives a very tiny weight fac- 
tor exp(— Mg++/Tc) 0.003 at Tc in comparison with 
unity factor of the vacuum state. That is to say, for the 
quenched glueballs, up to the critical temperature Tc, 
the contribution of higher spectral components beyond 
the vacuum to the spectral function are much smaller 
than the statistical errors (the relative statistical errors 
of the thermal glueball correlators are always a few per- 
cent) and can be neglected. As a result, the function 
form of p{i-u) in Eq. [15] can be a good approximation for 
the spectral function of glueballs at least up to T^. Ac- 
cordingly, considering the finite extension of the lattice in 
the temporal direction, the function form of the thermal 
correlators can be approximated as 



cosh(M„(l/(2r)-0) 
sinh(Af„/(2r)) ■ 



(17) 



which is surely the commonly used function form for the 
study of hadron masses at low temperatures on the lat- 
tice. As is always done, the glueball masses M„ derived 
by this function are called the pole masses in this work. 



A. Results of the single-cosh fit 

Even though the above discussion are based on the 
weak-interaction approximation for the hadronlike modes 
below Tc, we would like to apply Eq. [17] to analyzing 
the thermal correlators all over the temperature in con- 
cern. The interest of doing so is twofold. First, the ther- 
mal scattering of the glueball-like modes would result in 
a mass shift, say the deviation of the pole mass from 
the glueball mass at zero-temperature, which reflects the 
strength of the interaction at different temperature. Sec- 
ondly, the breakdown of this function form would signal 
the dominance of new degrees of freedom instead of the 
hadronlike modes in the thermal states. 

In practice, after the thermal correlators C{t, T) of the 
optimal operators are obtained according to the steps 
described in Sec. 11(B), the pole masses of the ground 
state (or the lowest spectral component) can be extracted 
straightforwardly. First, for each i?^'-^ channel and at 
each temperature T, the effective mass Mcs{t) as a func- 
tion of t is derived by solving the equation 



C{t + 1, T) _ cosh((t + 1 - Ntl2)atM,a{t)) 



C{t,T) 



cosh((i - Nt/2)atMcs{t)) 



(18) 



Secondly, the effective masses are plotted versus t and 
the plateaus give the fit windows [^1,^2]- Finally, the 
pole masses of the ground states are obtained by fitting 
C{t, T) through a single-cosh function form. As a con- 
vention in this work, we use Mq to represent the mass of 
a glueball state in the physical units and M to represent 
the dimensionless mass parameter in the data processing 
with the relation AI = Meat- 
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FIG. 4: Effective masses at different temperatures in A^^ 
channel. Data points are the effective masses with jackknife 
error bars. The vertical lines indicate the time window [ti, t2] 
over which the single-cosh fittings are carried out, while the 
horizontal lines illustrate the best-fit result of pole masses (in 
each panel the double horizontal lines represent the error band 
estimated by jackknife analysis) 



In Fig.[31 Fig.[Sl Fig.[Sl and Fig.[7]are shown the effec- 
tive masses with jackknife errors at various temperatures 
in Af^, ^r^j and channels, respectively. The 

vertical lines indicate the time window [ti, ^2] over which 
the single-cosh fittings are carried out, while the hori- 
zontal lines illustrate the best-fit result of pole masses 
(in each figure panel the double horizontal lines give the 
error band estimated by jackknife analysis). These fig- 
ures exhibit some common features: At the tempera- 
tures below Tc {Nt = 128, 80, 40), the effective mass 
plateaus show up almost from right the beginning of t, 
as it should be for the optimal glueball operators, while 
at T > Tc {Nt = 36, 24, and 20), the plateaus appear 
later and later in time, and even do not exist at Nt = 20 
(T — 1.90Tc). This observation can be interpreted as 
follows. Since the effective masses are calculated based 
on Eq. I17[ the very early appearance of the plateaus be- 
low Tc implies that the thermal correlators C{t, T) of the 
optimal operators are surely dominated by the ground 
state and can be well described by the function form of 
Eq. [T7] In other words, the picture of weakly interacting 
glueball-like modes makes sense for the state of matter 
below Tc. While at T > Tc, the later appearance and 
the narrower size plateaus signal that the picture of the 
state of matter is distinct from that at T < Tc. How- 
ever, because of the existence of effective mass plateaus 



7 



■ M,[f(N 



5 lU 15 20 25 30 



5 lU 15 20 25 30 



l\l^N,=l28) . 



5 10 15 20 25 30 



5 10 15 20 25 30 





'M,I^N,=40)' > 











2 4 6 8 10 12 14 16 18 



0.15 
0.1 
0.05 



M,I^N,=2 


1) ^ 















:36)' . » . 









2 4 6 8 10 12 14 16 18 



M^„<N,=20) . 







1,=40)' . . . ' 








U 2 4 6 8 10 12 14 16 18 



M,.ff(N|=2t) 





>N,=36) 











2 4 6 8 10 12 14 16 18 



M^„<N,=20) . 



2 4 6 8 10 12 



2 4 6 8 10 



2 4 6 8 10 12 



2 4 6 8 10 



FIG. 5: Similar to Fig. [H but in + channel. 



FIG. 7: Similar to Fig.gl but in T^'^ channel. 
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FIG. 6: Similar to Fig. H but in E++ channel. 



up to r - 1.58rc(iVt = 24), the possibility that glueball- 
like modes survive at this high temperature cannot be 
excluded. 

The pole masses in all 20 R^'-^ channels are extracted 
in units a^^ at all temperatures and are shown in Ta- 
ble HVl Specifically, with the lattice spacing determined 
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FIG. 8: The T-dependence of pole masses A~+ 
and glueballs. 
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in Sec. II, the pole masses of A^^ , £■++ and T. 

at T ~ in physical units are M^++=1.576(22)GeV, 
Af^-+=2.488(34)GeV, Me^^ ^ Mj,++=2.364(ll)GeV, 
respectively, which are in agreement with that of previ- 
ous studies |2l| - [28| . From the table, one can see that the 
behaviors of the pole masses with respect to the tempera- 
ture in all 20 channels are uniform: the pole masses keep 
almost constant with the temperature increasing from 
0.30Tc to right below (0.95Tc), and start to reduce 
gradually when T > Tc- When T increases up to l.OOTc, 
the pole masses cannot be extract reliably through the 
single-cosh fit for the lack of clear effective mass plateaus. 
Figure [5] iUustrates these behavior of pole masses in A^^, 
E++ and T++ channels. 
These results imply that glueballs can be very stable 
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TABLE IV: The pole 


xicLSses ( 


in units of ^) 


in all the 20 R^^ 


channels are extracted at all the temperatures. 


dPC 

R 


128 




80 




60 




48 


44 


40 




36 


32 


28 


24 




0.140( 


2) 


0.144( 


3) 


0.144( 


2) 


0.143( 3) 


0.140(2) 


0.140( 


3) 


0.132( 4) 


0.126( 2) 


0.122( 4) 


0.116( 3) 


1 _| 

At 


0.441( 


3) 


0.435( 


3) 


0.434( 


5) 


0.437( 4) 


0.432( 4) 


0.435( 


5) 


0.399( 6) 


0.322( 9) 


0.267(16) 


0.241(13) 




0.221( 


3) 


0.225( 


2) 


0.222( 


2) 


0.225( 2) 


0.218( 3) 


0.222( 


2) 


0.183( 5) 


0.174( 3) 


0.155( 4) 


0.146( 4) 




0.475( 


6) 


0.453( 


8) 


0.447( 


9) 


0.464( 7) 


0.473( 6) 


0.468( 


6) 


0.426(12) 


0.417(10) 


0.287(19) 


0.253(18) 


At+ 


0.323( 


4) 


0.327( 


4) 


0.326( 


4) 


0.330( 2) 


0.326( 4) 


0.332( 


3) 


0.282( 7) 


0.249( 8) 


0.224( 9) 


0.208( 9) 


.J 

At 


0.302( 


5) 


0.308( 


3) 


0.308( 


5) 


0.312( 3) 


0.312( 5) 


0.308( 


6) 


0.268( 6) 


0.241( 7) 


0.220( 8) 


0.201( 6) 


AI+ 

1 


0.450( 


5) 


0.449( 


7) 


0.446( 


5) 


0.440( 6) 


0.452( 4) 


0.448( 


5) 


0.396(10) 


0.340(11) 


0.330(12) 


0.250(14) 


Ar 


0.387( 


3) 


0.388( 


3) 


0.385( 


4) 


0.390( 5) 


0.376( 4) 


0.375( 


4) 


0.354( 7) 


0.293( 7) 


0.268(10) 


0.214( 9) 


£;++ 


0.210( 


1) 


0.205( 


1) 


0.207( 


1) 


0.209(2) 


0.206(1) 


0.189( 


4) 


0.167( 4) 


0.153( 3) 


0.143( 3) 


0.139( 2) 




0.401( 


2) 


0.403( 


2) 


0.401( 


2) 


0.394( 4) 


0.400( 2) 


0.395( 


3) 


0.375( 4) 


0.311( 6) 


0.261( 7) 


0.230( 7) 


E-+ 


0.273( 


1) 


0.266( 


1) 


0.264( 


2) 


0.273( 2) 


0.275( 1) 


0.262( 


2) 


0.218( 4) 


0.196( 4) 


0.183( 4) 


0.181( 4) 


E-- 


0.374( 


1) 


0.368( 


2) 


0.360( 


2) 


0.361( 3) 


0.363( 3) 


0.352( 


4) 


0.308( 8) 


0.262( 6) 


0.231( 6) 


0.213( 6) 


T++ 


0.327( 


2) 


0.326( 


4) 


0.327( 


2) 


0.334( 2) 


0.331( 2) 


0.312( 


5) 


0.287( 7) 


0.266( 3) 


0.227( 6) 


0.215( 4) 


r+- 


0.278( 


1) 


0.274( 


2) 


0.265( 


3) 


0.278( 2) 


0.281( 1) 


0.261( 


3) 


0.207( 6) 


0.199( 2) 


0.181( 4) 


0.175( 2) 




0.372( 


2) 


0.377( 


4) 


0.371( 


3) 


0.380( 2) 


0.374( 2) 


0.370( 


3) 


0.331( 5) 


0.289( 7) 


0.248( 7) 


0.230( 5) 




0.350( 


4) 
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2) 
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0.205( 
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1) 
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3) 


0.160( 3) 
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2) 


0.317( 


2) 
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4) 
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0.303( 


5) 
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2) 


0.260( 


3) 


0.264( 


2) 
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2) 


0.240( 3) 
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0.187( 4) 
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2) 


0.358( 


3) 
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3) 


0.358( 4) 


0.367( 2) 


0.353( 


5) 


0.282(13) 


0.254( 6) 


0.235( 6) 
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below Tc and survive up to l.GTc- This coincides with 
the thermal properties of heavy quarkonia observed by 
model calculation and lattice numerical studies [^-[S lla - 
[TqI . liol l4l| , but different from the observation of a pre- 
vious lattice study on glueballs where the observed pole- 
mass reduction start even at T ~ O.ST^ ISal. 



B. Breit-Wigner analysis 

In the single-cosh analysis, it is seen that, when the 
temperature increases up to Tc, the thermal correlators 
can be well described by Eq. [17] and the pole masses of 
glueballs are insensitive to T. This is in agreement with 
the picture that the state of matter below Tc are made up 
of weakly interacting glueball-like modes. When T > 
the thermal correlators deviate from Eq. [T7] more and 
more. This observation implies that the degrees of free- 
dom are very different from that when T < Tc- Theo- 
retically in the deconfined phase, gluons can be liberated 
from hadrons. However, the study of the equation of 
state shows that the state of the matter right above Tc 
is far from a perturbative gluon gas. In other words, the 
gluons in the intermediate temperature above Tc may 
interact strongly with each other and glueball-like res- 
onances can possibly be formed. Thus different from 
bound states at low temperature, thermal glueballs can 
acquire thermal width due to the thermal scattering be- 
tween strongly interacting gluons and the magnitudes of 
the thermal widths can signal the strength of these types 



of interactions at different temperatures. 

In order to take the thermal width into consideration, 
we also adopt the Breit-Wigner Ansatz, which is sug- 
gested by the pioneering work Ref. [s^, to analyze the 
thermal correlators once more. First, we treat thermal 
glueballs as resonance objects which correspond to the 
poles (denoted by = luq — iT) of the retarded and 
advanced Green functions in the complex oj.plane (note 
that conventionally in particle physics, a resonance pole 
is always denoted as M — iT /2 where M is the mass of 
the resonance and F is its width.) cjq is called the mass 
of the resonance glueball and F its thermal width in this 
work. Secondly, we assume that the spectral function 
p{uj) is dominated by these resonance glueballs. Thus 
the spectral function is parametrized as 



p{uj) = A((5r(w - Wo) — 5t{oj + uoq) + 
where 5, is the Lorentzian function 



5t{x) = — Im 

TT 



1 



1 



(19) 



(20) 



and " . . ." represents the terms of excited states. With 
this spectral function, the thermal glueball correlator 
G{t, T) can be expressed as 



C{t,T) 



duj cosh{uj{^ ~ t)) 
2^ 2sinh(^) 

— OO 

X 27r^ {5r{uj - ujq) - Sr{uj + ojq)) + ... .(21) 
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FIG. 9: Determinations of the fit range [ti,t2] in channel at Nt = 128, 36, and 20. In each row, ui^o°^\t) and F^^^^t) 

obtained by solving Eq. [23] are plotted by data points with jackknife error bars. [ti,t2] are chosen to include the time slices 
between the two vertical lines, where LOg^^^t) and r'"^^(t) show up plateaus simultaneously. The best-fit results of loo and F 
through the function gr{t) are illustrated by the horizontal lines. 



The integral on the right hand side of above equation, denoted by gr{t), can be calculated explicitly as 
I 



grit) = A 



Re 



cosh((a;o + jr)(2^ - t)) 



sinh( 



2T 



oo 

2ujqT ^ cos {2nnTt) 



(27rnT + r)^ 



— (n — > — n) 



(22) 



which can be used as the fit function to extract ujq and F 
from the thermal correlators obtained from the numerical 
calculation. Practically, the infinite series in the above 
equation is truncated by setting the upper limit of the 
summation to be 50, which is tested to be enough for all 
the cases considered in this work. 

In the present study, we carry out the Breit-Wigner 
analysis in A^'^ , ^ -E'^^j a-nd T^~^ channels, whose 
continuum correspondences are ^, and 2++. Al- 
though the variational method is exploited to enhance the 
contribution of the ground state to the thermal correla- 
tors, the contributions from higher spectral components 
cannot be eliminated completely. Therefore, the fit range 
must be chosen properly where the contribution of the 
ground state dominates. We take the strategy advocated 
in Ref.fssj as follows. For a given correlator C{t,T), the 
effective peak position LL)j^^\t) and the effective width 
p(cff)('^^ are obtained by solving the equations 

grjt) ^ C{t,T) 
grit + I) C(t + l,r)' 



gr(^ + l) ^ C{t + l,T) 

grit + 2) Cit + 2,T)' ^ ' 

The statistical errors oiJ^^\t) and F(°ff)(i) can be esti- 
mated through the jackknife analysis. Thus the fit range, 
denoted by [ti , ^2] , is chosen to be the time range where 

u]'^^\t) and T^'^^\t) show up plateaus simultaneously. 
For example, the procedure in T^"*" channel is illustrated 
in Fig. [S] for Nt = 128,36,20 (corresponding to the tem- 
perature T/Tc = 0.30,1.05,1.90), where the fit ranges 
[ti,t2\ are determined to include the time slices between 
the two vertical lines in each figure. 

After the fit ranges for all the thermal correlators are 
chosen, the jackknife analysis can be carried out straight- 
forward and the detailed procedures are omitted here. 
TableEl ED EU and Einishow the fit windows [tiM]. 
the chi-square per degree of freedom ^ jd.o.j , and the 
best-fit results of wq and F at various temperature in 
, A'{^ , E++, and T^^ channels. In almost aU the 
cases, the fit ranges start from ti = 1, 2, or 3, and last for 
quite a few time slices. This reflects that, as is expected, 
the optimal glueball operators couple almost exclusively 



10 



A A A&A 


* +'+ A ' 

A. 1 — A — 1 
T/Tp=1.0 - 

A A A A A - 


AAA AA^ 


T/Tj=I.O 

A A A 


\ 

AAA AAa 


E^^^^ 

T/Tj=1.0 


1 1 

AAA aAa 


' T*'* . A ,' 

T/Tj=i.O 

A A A A - 



0.5 1 1.5 2 0.5 1 1.5 2 
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TABLE V: The best-fit ujq and F of A++ channel at different 
T through the Breit-Wigner fit. Also listed are the fit window 
[ti,f2] and the chi-square per degree of freedom, jd.o.f . 



Nt 


T/Tc 




F 




x'/d.o.f 


128 


0.30 


0.142(2) 


0.008(10) 


(3, 6) 


1.266 


80 


0.47 


0.146(2) 


0.013(9) 


(2, 4) 


0.921 


60 


0.63 


0.144(2) 


0.008(3) 


(2, 7) 


0.135 


48 


0.79 


0.143(2) 


0.014(6) 


(2, 4) 


0.639 


44 


0.86 


0.142(2) 


0.004(3) 


(1, 7) 


0.758 


40 


0.95 


0.143(2) 


0.003(4) 


(4, 7) 


0.850 


36 


1.05 


0.146(2) 


0.028(4) 


(3, 6) 


0.960 


32 


1.19 


0.141(2) 


0.053(2) 


(1,4) 


0.393 


28 


1.36 


0.142(2) 


0.056(4) 


(2, 5) 


1.253 


24 


1.58 


0.143(1) 


0.059(3) 


(1, 4) 


0.302 


20 


1.90 


0.146(2) 


0.077(5) 


(2, 4) 


0.918 
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FIG. 11: F's are plotted versus T/Tc for A-+ E++ , 

and T^^ channels. The vertical lines indicate the critical 
temperature. 



TABLE VI: The best-fit a;o and T of A~+ channel at different 
T through the Breit-Wigner fit. Also listed are the fit window 
[ti,t2] and the chi-square per degree of freedom, /d.o.f. 



Nt 


T/T, 




F 




xVd.o.f 


128 


0.30 


0.226(2) 


0.006(5) 


(3, 9) 


0.509 


80 


0.47 


0.228(2) 


0.008(4) 


(2, 6) 


0.640 


60 


0.63 


0.227(1) 


0.013(5) 


(2, 7) 


0.216 


48 


0.79 


0.228(2) 


0.012(5) 


(2, 6) 


0.177 


44 


0.86 


0.229(2) 


0.013(4) 


(2, 8) 


0.184 


40 


0.95 


0.224(2) 


0.004(6) 


(3, 6) 


0.549 


36 


1.05 


0.221(2) 


0.037(3) 


(1,8) 


0.935 


32 


1.19 


0.219(2) 


0.047(3) 


(1,6) 


0.250 


28 


1.36 


0.211(3) 


0.068(6) 


(2, 5) 


0.091 


24 


1.58 


0.208(3) 


0.089(7) 


(2, 4) 


0.003 


20 


1.90 


0.211(3) 


0.099(9) 


(2, 6) 


0.083 



to the lowest spectral components after the implemen- 
tation of the variational method. All the /d.o.f^s are 
^ 0(1) or even smaller, which reflect the reliability of 
the fits. 

The main features of the best fit wq and F based on 
Breit-Wigner Ansatz are described as follows: 

• The peak positions ujq of the spectral functions p{uj) 
are insensitive to the temperature in all the consid- 
ered channels. In particular, the wq in A'^'^ chan- 
nel keeps almost constant all over the temperature 
range from O.SOTc to 1.90Tc. In the other three 
channels, the wo's do not change within errors be- 
low Tc, but reduce mildly with the increasing tem- 
perature above Tc- The reduction of wq at the high- 
est temperature T — 1.90Tc is less than 5% in these 
three channels. 

• In all four channels, the thermal widths F are small 
and do not vary much below Tc, but grow rapidly 
with the increasing temperature when T > Tc. Be- 
low Tc, the thermal widths are of order F ~ 5% or 
even smaller (especially for the Af~^ F is consistent 
with zero). The thermal widths increase abruptly 



when the temperature passes Tc and reach values 
- cjo/2 at T = l.QOTc. 

These features can be seen easily in Fig. [TU] and [TTl 
where the behaviors of ujq and F with respect to the 
temperature T are plotted for all four channels. The 
line_shapes of the spectral functions with the best-fit pa- 
rameters at different T are shown in Fig. [12] for 
E'^'^, and Tj^"*" channels (we do not plot the spec- 
tral function of A^^ channel due to the small thermal 
widths). 



IV. SUMMARY AND DISCUSSIONS 

On 24"^ X Nt anisotropic lattices with the anisotropy 
^ = 5 at the gauge coupling /3 = 3.2, the thermal 
glueball correlators are calculated in a large tempera- 
ture range from O.SOTc to 1.90Tc, which are realized by 
varying Nt to represent different temperatures. Based 
on the lattice spacing Qs = 0.0878(4) fm determined by 
— (410(20) MeV), the spatial extension of the lattices 
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FIG. 12: Plotted are the spectral function p{ui) at T /Tc — 0.30, 0.95,1.05 and 1.58 with the best-fit parameters. Panel (a)_(c) 
are for A^^ , and T^^ channels, respectively. 



TABLE VII: The best-fit ljq and F of E++ channel at different 
T through the Breit-Wigner fit. Also listed are the fit window 
[ii,i2] and the chi-squarc per degree of freedom, /d.o.f. 



Nt 




[Jo 


F 




xVdof 


128 


0.30 


0.212(1) 


0.012(4) 


(2, 5) 


0.274 


80 


0.47 


0.211(1) 


0.006(3) 


(2, 8) 


0.616 


60 


0.63 


0.212(1) 


0.010(3) 


(2, 9) 


0.844 


48 


0.79 


0.213(1) 


0.011(3) 


(2, 5) 


0.206 


44 


0.86 


0.211(1) 


0.011(3) 


(2, 8) 


1.268 


40 


0.95 


0.207(1) 


0.022(3) 


(2, 6) 


0.250 


36 


1.05 


0.205(2) 


0.034(2) 


(1,8) 


0.183 


32 


1.19 


0.200(1) 


0.049(2) 


(1,6) 


0.478 


28 


1.36 


0.191(2) 


0.067(4) 


(2, 6) 


0.297 


24 


1.58 


0.189(2) 


0.083(5) 


(2, 4) 


0.253 


20 


1.90 


0.196(2) 


0.091(5) 


(2, 4) 


0.046 



TABLE VIII: The best-fit wo and F of T++ channel at dif- 
ferent T through the Breit-Wigner fit. Also listed are the fit 
window [ti,t2] and x^/d-o.f. 



Nt 






F 






128 


0.30 


0.210(1) 


0.008(2) 


(3,10) 


0.442 


80 


0.47 


0.213(1) 


0.009(3) 


(3, 9) 


0.696 


60 


0.63 


0.213(1) 


0.012(3) 


(3, 7) 


0.326 


48 


0.79 


0.210(1) 


0.007(4) 


(3, 6) 


0.288 


44 


0.86 


0.214(1) 


0.012(2) 


(2, 6) 


0.437 


40 


0.95 


0.208(1) 


0.023(1) 


(1,7) 


0.743 


36 


1.05 


0.204(1) 


0.039(2) 


(1,7) 


0.606 


32 


1.19 


0.199(1) 


0.047(1) 


(1,6) 


0.527 


28 


1.36 


0.196(2) 


0.064(3) 


(2, 5) 


0.022 


24 


1.58 


0.194(1) 


0.077(4) 


(2, 7) 


0.119 


20 


1.90 


0.196(2) 


0.093(4) 


(2, 5) 


0.027 



are estimated to be (2.1 fm)'^, which is large enough to 
be free of the finite volume effects. On the other hand, 
because of the large anisotropy, there are enough data 
point in the temporal direction for the thermal correla- 
tors to be analyzed comfortably even at the highest tem- 
perature T ~ 2Tc concerned in this work. With the im- 
plementation of the smearing scheme and the variational 



method, we can construct the optimal glueball operators 
in all the symmetry channel, which couple mostly to the 
lowest-lying states (or more precisely, the lowest-lying 
spectral components). As a result, the thermal corre- 
lators of these operators can be considered to be con- 
tributed dominantly from these lowest-lying states. The 
thermal correlators are analyzed based on two ansatz, 
say, the single-cosh function form and the Breit-Wigner 
Ansatz. In Table IIXI and Table |Xl the "pole masses" 
Mq obtained by single-cosh analysis and (wq, F) obtained 
based on the Breit-Wigner Ansatz are combined together 
for comparison (all the data are converted into physical 
units). 

The most striking observation from the single-cosh 
analysis is that, in all 20 R^^ channels, the best-fit pole- 
masses Mq are almost constant within errors from the 
low temperature up to right below the critical tempera- 
ture Tc- This is what should be from the point of view of 
deconfinement phase transition of QCD: Since below Tc 
the system is in the confinement phase, the fundamen- 
tal degrees of freedom must be hadrons. Above Tc, the 
reduction of the pole masses does signal the QCD transi- 
tion, after which the state of the matter is very different 
from that below Tc- However, the existence of effective 
mass plateaus, from which the pole masses are extracted, 
also implies that color singlet objects, the glueball-like 
modes, can also survive at the intermediate temperature 
above Tc- The results of the Breit-Wigner fit are consis- 
tent with this picture. In the Breit-Wigner Ansatz, ther- 
mal widths F are introduced to glueball states to account 
for the effects of finite temperature, such as the thermal 
scattering and the thermal fluctuations. As shown in Ta- 
ble IIXI and |Xl below Tc (or in the confinement phase) , 
the best-fit wq's are very close to the pole masses, and 
the thermal widths F are very tiny and are always of a 
few percent of wq. This means the glueball states are 
surely stable in the confinement phase and the thermal 
interaction among them are weak. With the temperature 
increasing above Tc, while the temperature dependence of 
wo's is very mild, the thermal widths F grow rapidly and 
reach values of roughly half of cj's at T ^ 1.9Tc. This 
clearly reflects that glueballs act as resonances are un- 
stable more and more, and the reduction of pole masses 
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TABLE IX: The "pole masses" Ma obtained by single-cosh analysis and (tJo,r) obtained based on the Breit-Wigner ansatz 
are combined together for comparison. Listed in the table are the results in A^^ and A^^ channels (all the data are converted 
into the physical units). 





AT* 


Nt 




mcpeV] 


a;o[GeV] 


r[GeV] 


mG[GeV] 


a;o[GeV] 


r[GeV] 


128 


0.30 


1.576(22) 


1.602(14) 


0.091(113) 


2.488(31) 


2.549(17) 


0.069(52) 


80 


0.47 


1.621(29) 


1.644(22) 


0.145(100) 


2.533(24) 


2.560(17) 


0.086(44) 


60 


0.63 


1.627(18) 


1.621(16) 


0.098(38) 


2.499(27) 


2.559(16) 


0.147(54) 


48 


0.79 


1.616(21) 


1.612(26) 


0.156(67) 


2.533(23) 


2.564(20) 


0.135(55) 


44 


0.86 


1.577(25) 


1.598(17) 


0.045(34) 


2.454(34) 


2.577(18) 


0.144(48) 


40 


0.95 


1.576(39) 


1.621(20) 


0.034(46) 


2.499(25) 


2.525(26) 


0.042(71) 


36 


1.05 


1.486(43) 


1.638(28) 


0.315(48) 


2.060(48) 


2.490(23) 


0.413(32) 


32 


1.19 


1.418(21) 


1.588(25) 


0.586(28) 


1.959(37) 


2.464(20) 


0.529(28) 


28 


1.36 


1.373(48) 


1.599(26) 


0.619(43) 


1.745(43) 


2.375(28) 


0.768(71) 


24 


1.58 


1.306(41) 


1.613(28) 


0.664(37) 


1.644(45) 


2.308(32) 


0.999(83) 


20 


1.90 




1.642(32) 


0.873(59) 




2.380(35) 


1.114(99) 



TABLE X: The pole masses Mg obtained by single-cosh analysis and (a;o,r) obtained based on the Breit-Wigner Ansatz are 
combined together for comparison. Listed in the table are the results in E'^^ and T^"*" channels (all the data are converted 
into physical units). 



E++ 


rn+ + 


Nt 


T/T, 


mG[GeV] 


a;o[GeV] 


r[GeV] 


mG[GeV] 


ajo[GeV] 


r[GeV] 


128 


0.30 


2.364(11) 


2.385(12) 


0.140(42) 


2.308(14) 


2.363(10) 


0.091(25) 


80 


0.47 


2.308(13) 


2.368(12) 


0.069(29) 


2.353(15) 


2.387(10) 


0.105(32) 


60 


0.63 


2.330(11) 


2.383(12) 


0.116(32) 


2.319(13) 


2.396(10) 


0.140(32) 


48 


0.79 


2.353(19) 


2.393(15) 


0.129(34) 


2.308(15) 


2.362(14) 


0.083(43) 


44 


0.86 


2.319(15) 


2.379(12) 


0.119(32) 


2.330(21) 


2.405(10) 


0.136(26) 


40 


0.95 


2.263(14) 


2.327(15) 


0.247(38) 


2.218(16) 


2.344(11) 


0.259(16) 


36 


1.05 


1.880(41) 


2.305(17) 


0.382(23) 


1.925(33) 


2.298(14) 


0.437(17) 


32 


1.19 


1.722(35) 


2.247(16) 


0.549(20) 


1.801(25) 


2.244(11) 


0.532(15) 


28 


1.36 


1.610(31) 


2.155(19) 


0.754(43) 


1.666(23) 


2.205(17) 


0.717(36) 


24 


1.58 


1.565(23) 


2.132(21) 


0.937(55) 


1.610(27) 


2.184(16) 


0.870(41) 


20 


1.90 




2.201(23) 


1.023(60) 




2.209(20) 


0.935(48) 



above Tc can be taken as the effect of these growing ther- 
mal widths. 

To summarize, in pure gauge theory, the state of mat- 
ter is dominated by weakly interacting hadronlike states 
below Tc] when T > T^, glueball states survive as res- 
onancelike modes up to a temperature T ~ 1.9Tc with 
their thermal widths growing with increasing T, which 
implies that in this intermediate temperature range, glue- 
balls are unstable and may decay into gluons, and re- 
versely gluons also interact strongly enough to form 
glueball-like resonances. The two procedure may reach 
the thermal equilibrium at a given temperature, such 
that the gluon degree of freedom become more and more 
important with T increasing. At very high temperature, 
the glueball-like resonances may disappear finally and the 
state of matter can thereby be described by a perturba- 
tive gluon plasma. This picture is coincident with the ob- 
servations both in the study of equation of state of QCD 
and the thermal properties of heavy quarkonia. On the 
other hand, the surprising results of RHIC experiments 
mav also suDDort this oicture to some extent. First, the 



data of RHIC experiments are well described by the hy- 
drodynamical model[42,]. Secondly, the investigation of 
elliptic flow data using a Boltzmann-type equation for 
gluon scattering is not consistent with the perturbative 
QCD apparently [i^. So the quark-gluon plasma at the 
RHIC temperature is most likely a strongly interacting 
system. 
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